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Abstract 



^ ' The scale of mass generation for fermions (including neutrinos) and the scale for electroweak 

' , symmetry breaking (EWSB) can be bounded from above by the unitarity of scattering involving 

longitudinal weak gauge bosons or their corresponding would-be Goldstone bosons. Including the 
(~| , exact n-body phase space we analyze the 2 — > n (n ^ 2) processes for the fermion-(anti)fermion 

scattering into multiple gauge boson final states. Contrary to naive energy power counting, 
we demonstrate that as n becomes large, the competition between an increasing energy factor 
' and a phase-space suppression leads to a strong new upper bound on the scale of fermion mass 

generation at a finite value n = ng, which is independent of the EWSB scale, v — (V2Gi?)~2. 
For quarks, leptons and Majorana neutrinos, the strongest 2 — > n limits range from about 3 TeV 
to 130— 170 TeV (with 2 < < 24), depending on the measured fermion masses. Strikingly, 
given the tiny neutrino masses as constrained by the neutrino oscillations, neutrinoless double- 
beta decays and astrophysical observations, the unitarity violation of I'j^V]^ nW^ scattering 
actually occurs at a scale no higher than ~ 170 TeV. Implications for various mechanisms of 
neutrino mass generation are analyzed. On the other hand, for the 2 n pure Goldstone- 
boson scattering, we find that the decreasing phase space factor always dominates over the 
growing overall energy factor when n becomes large, so that the best unitarity bound on the 
scale of EWSB remains at n ~ 2 . 
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1. Introduction 



Despite the astonishing success of the standard model [1, 2] for electroweak and strong interactions 
over the past thirty years, the origin and scales of mass generation remain, perhaps, the greatest 
mystery. The standard model (SM) contains nineteen free input parameters in total, among which 
twelve can be expressed in terms of unpredicted mass-eigenvalues (two weak gauge boson masses, 

one Higgs mass, six quark masses and three lepton masses) with one good equality^ ^ ~ 1 

v/V2 

IJle _R 

and a huge hierarchy ~ 3 x 10 as fixed by experimental data. There are four additional 

mt 

inputs (three mixing angles and one CP-phase) in the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
[3] coming from diagonalizing the quark mass matrices, so that sixteen out of the nineteen free 
parameters^ arc due to our lack of knowledge about the origin of mass generation. The recent 
exciting neutrino oscillation experiments [4, 5, 6, 7, 8] further point to three extra neutrino masses^ 

at the scale of 0.05 eV, giving another huge hierarchy: — ~ 3 X 10"^^. In the SM, a single 

mt 

hypothetical Higgs boson [18] is assumed to generate masses for all gauge bosons, quarks and leptons, 
which, however, leaves all mass scales and their hierarchies fixed by empirical inputs, without real 
explanation. Furthermore, such a Higgs boson cannot generate neutrino masses without losing 
renormalizability (given the observed SM particle spectrum) [19] or adding extra new particles (to 
retain the renormalizability) [20, 21]. Various interesting extensions of the SM have been sought, 
ranging from compositeness [22, 23, 24], low or high scale supersymmetry [25, 26], grand unification 
[27], to extra dimensions and deconstructions [28, 29], and string/M theories [30]. Each of them 
overcomes some unsatisfactory aspects of the SM while adding many new unknowns. Among so 
many kinds of possible new physics, there is no doubt that understanding the mass generation 
holds the greatest promise to unravel immediate new clues beyond the SM. Fortunately, on-going 
and upcoming high energy experiments will actively explore the scales of mass generation (and 
the associated mixings) for weak gauge bosons, quarks, leptons and neutrinos, either directly or 
indirectly, in this decade. 

While we still seem a long way from finding an imagined theory of everything with only a single 

^Here v = (v^Gf) ^ — 246 GeV is determined by the Fermi constant Gf- 

^The remaining three parameters are two gauge couplings {os, aem) and one strong CP phase. 

^For Majorana neutrinos there are three mixing angles [9] and three CP phases in the leptonic sector [10]. The 
global analysis [11] of all the current oscillation data strongly disfavors an extra singlet sterile neutrino. The upcoming 
Fermilab MiniBooNE experiment [12] will resolve the controversy with LSND [13]. 

*The recent astrophysics analyses of the WMAP (Wilkinson Microwave Anisotropy Probe) data [14], in conjunc- 
tion with the 2dF Galaxy Redshift Survey [15], give an impressive 95%C.L. upper bound, ^jTUvj < 1.01 eV [16]. 
Combining the WMAP, SDSS (Sloan Digital Sky Survey) and Lyman-a Forest data results in a better limit of 
'^jTu^j ^ 0.65 eV (95%C.L.) [17]. 
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input parameter to explain the world, we wish to ask a much more modest "bottom-up" question: 
How do wc use high energy scattering processes as a fundemental means to probe the scales of 
mass generation in a model-independent way? What is the energy scale at which the new physics 
underlying the relevant mass generation has to show up? In particular, fermion mass generation 
can be independent of the electrowcak symmetry breaking (EWSB) that generates the gauge boson 
masses of the and Z^; since the fermion and gauge boson masses need not originate from the 
vacuum expectation value (VEV) of the same Higgs doublet (as assumed in the SM)^, we would 
hope to probe the scales for fermion mass generation without any ad hoc theory assumption. The 
SM gauge symmetry forbids all bare mass terms for gauge bosons and fermions. If any observed 
particle mass is put in by hand in the SM Lagrangian, the theory is no longer renormalizable and we 
expect the high energy scattering amplitudes involving this particle to grow with energy and violate 
unitarity at a certain scale above which the new physics has to show up. In this way the unitarity 
of the S'-matrix puts a universal upper bound on the scale of the corresponding mass generation. 

In the standard model EWSB sector for generating the W/Z masses, once the neutral physical 
Higgs boson is removed, the 2 — > 2 quasi-clastic scattering amplitudes of longitudinally polarized 
W/Z's exhibit the worst high energy behavior [34, 35, 36], proportional to the square of the cm. 
energy -y/s , leading to partial wave unitarity violation at a scale 



VSttv ~ 1.2 TeV, (1.1) 



similar to the classic analyses of the upper bound on the mass of the Higgs boson [35, 37, 38] where 
tree-unitarity is violated at a critical Higgs mass M|^'* ~ ~ 1 TeV ^. Eq. (1.1) puts a model- 

independent upper limit on the scale of EWSB and justifies TeV energy scale for the Large Hadron 
Collider (LHC) at CERN. Appelquist and Chanowitz derived an analogous upper bound on the 
scale of fermion mass generation by considering the 2^2 inelastic channels of fermion-antifermion 
scattering /±/± ^ Vp vp (V = W^,Z<>) [43], 

. (1.2) 

where Nc = 3(1) for quarks (leptons). This bound indicates that the upper limit on the scale of 
fermion mass generation is proportional to the inverse power of the fermion mass itself and is thus 

^For instance, in the Top-color models [31, 22], the top-condensate Higgs (which gives the observed large top- 
quark mass) only changes the EWSB vacuum expectation value {v ~ 246 GeV) by a few percent and thus has 
little to do with the EWSB. In the traditional dynamical models, the EWSB is caused by technicolor (TC) gauge 
interactions [32, 22], while the fermion masses are generated by technifermion condensates via additional extended 
technicolor (ETC) [33, 22], independent of EWSB. 

^Subsequent studies and applications may be found, for instance, in Refs. [39, 40, 41, 42] and references therein. 
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an independent bound from the customary unitarity bound on the EWSB scale in Eq. (1.1). Also, 
for all the SM fermions (except the top quark), the bound (1.2) is much weaker than Eq. (1.1). 
However, this traditional wisdom was recently challenged by an interesting study [44], where it was 
noted that for inelastic scattering with a multi-gauge-boson final state, // nV^ (n > 2), the 
n-body phase space integration contributes a nontrivial energy factor s"^^ to further enhance the 
cross section (in addition to the energy dependence of from the squared amplitude), so that the 
unitarity bound goes like 

— > V, (for large), (1.3) 

ijifj 

which could be pushed arbitrarily close to the weak scale v and thus become independent of the 
fermion mass m/ for large enough n. Prom this, Ref. [44] deduced that the fermion-antifermion 
scattering into weak bosons does not reveal an independent new scale for fermion mass generation. 
This is very counter-intuititive since the following kinematic condition must hold for any 2 — ^ n 
scattering with n weak gauge bosons in the final state, 

^s > nMw(z) - 3^ ' =^ V ^ 3 ■ ^^'^^ 

This shows that as n becomes large, the required cm. energy ^/s has to be arbitrarily above the 
weak scale v, so that the unitarity violation scale, y/s = E*, has to grow at least linearly with n, 
instead of decreasing to v , contrary to the power counting argument in Eq. (1.3). 

The observation in (1.4) is completely general, valid for any 2 — > n process including both // 
nV^ and V^^V^^ nVf^. Thus, a similar contradiction may possibly occur for the VlVl uVl 
scattering. We note that the n-body phase space is the same as before but the squared amplitude 
now behaves as s^/u^" , so by naive power-counting the inelastic unitarity bound goes as, 

E'^ ~ v{Cq)^ — y V, (for large), (1.5) 

where Cq > 1 is a dimensionless constant determined by the phase-space and the square of the 
S-matrix element. For n = 2, we know from Eq. (1.1) that Cq 2± (\/87r)^ — 632. Hence, as n 
becomes large, we would naively expect E^^ to approach v ~ 246 GeV, significantly below the 
traditional 2^2 bound E^^ ~ 1.2 TeV in Eq. (1.1). If this argument holds, we would deduce the 
new physics scale associated with EWSB to be less than ~ 250 GeV, which could manifest itself as 
a light Higgs boson or something similar and might already be partly accessible at the on-going 
Tevatron Run-2. 
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It is the above general observation based upon the kinematic condition (1.4) that leads us to 

believe that the energy power counting arguments in Eqs. (1.3) and (1.5) cannot be true and a new 
resolution must be sought. Therefore, in this study, we systematically and quantitatively analyze 
the 2 ^ n inelastic channels, with n > 2, for both // nVi and VlVl nVi scattering. 
Our findings show that the exact n-body phase space contributes a non-trivial dimensionless factor 
which sufficiently suppresses the effect of the energy growth given by Eqs. (1.3) and (1.5). As 
a consequence, the kinematic condition (1.4) is satisfied for any large n value, and the scale for 
fermion mass-generation is proven to be independent of the EWSB scale. 

Section 2 establishes the necessary concepts and formalism, and reviews the customary unitarity 
bounds from 2^2 scattering. We first define the scale of mass generation in connection with the 
unitarity violation scale. We then discuss the equivalence theorem, give a general power counting 
formula and apply it to count the energy dependence of the relevant scattering amplitudes for our 
later analyses. We present a systematical derivation of the general unitarity condition for 2 ^ n 
scattering, extending that of Ref. [44]. We then summarize the customary unitarity limits from 
2 — > 2 scattering for the scale of EWSB and the scales of mass generation for Dirac fermions and 
Majorana neutrinos. 

In Sec. 3 we return to the new puzzle presented around Eqs. (1.3)-(1.4) and resolve it by using an 
exact expression for the n-body phase space. With estimated amplitudes for the 2 — > n scattering, 
we analyze the unitarity violation scales for the scatterings ^1^2 — n where the ^ denotes 
quarks, leptons, Majorana neutrinos, or longitudinal weak bosons. 

In Sec. 4, 5, and 6 we further quantify the estimates in Sec. 3 for quarks and leptons, for Majorana 
neutrinos, and for EWSB, respectively. In Sec. 4 we first compute the scattering amplitude for 
// — nVl (mr^) from the universal contact interaction. This gives slightly weakened limits for the 
unitarity violation, i.e., the estimates of Sec. 3 are relaxed by a factor of 2 or so. We then examine 
the effect of including model-dependent EWSB contributions and note that the effects decrease with 
the increase of n, and, in particular, we explicitly show that for a SM (or SM-Uke) EWSB sector 
the effect on 2 ^ 3 scattering is to enhance the bounds slightly. Thus, we find that for all light 
SM fermions (except possibly the top quark) the 2 — n unitarity limits on the scale of fermion 
mass generation become the strongest at a value n = Ug > 2 and are substantially tighter than the 
classic Appelquist-Chanowitz bound [43]. 

In Sec. 5, we present a quantitative calculation of the scattering amplitude for Majorana neu- 
trinos i^i^i']^ nVl . This relaxes the estimate in Sec. 3 by a factor of about 1.3 — 1.6. Thus 
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we demonstrate that, without assuming' a priori that the same SM Higgs doublet used for EWSB 

participates in neutrino mass generation, the unitarity of the scattering i^iVj^ nV^ puts an upper 
bound on the Majorana neutrino mass generation at a scale substantially below M^yrp ,^ around 
130 — 170 TcV, with ~ 20 — 24 . We also discuss the implications of such low-scale unitarity 
violation for various mechanisms of Majorana neutrino mass generation. 

Finally, in Sec. 6, we explicitly analyze the scattering V^^Vj^^ nVf^ (vr^ivr^^ nvr"), with 
n ^ 4 , and show that as n increases the upper bound on the EWSB scale is not driven towards the 
scale V ~ 246 GeV [cf. Eq. (1.5)], rather, it becomes weaker so that the best constraint still occurs 
at Us = 2, in agreement with Sec. 3. Including our earlier general estimates from Sec. 3 for large n 
values, we conjecture that this feature holds for any value of n ^ 2 . 

We conclude in Sec. 7 and derive the expression for n-body phase space in Appendix A. For 
completeness, in AppendixB we also refine the unitarity conditions of Sec. 2.3.1. 



2. Scales of Mass Generation vs. Unitarity Bounds 

In this section, we start by introducing necessary concepts and formalism for the gauge and fermion 
mass terms and for the 2 ^ n (n > 2) unitarity limit. Then we briefly review the traditional 2 — > 2 
unitarity bounds derived for EWSB and fermion (neutrino) mass generation. 

2.1. Defining tiie Scale for Mass Generation 

It turns out that the SM Lagrangian, with the physical Higgs boson removed (often called the 
Higgsless SM) and with all observed particle masses put in by hand, exhibits a nonlinearly realized 
gauge symmetry [47, 48, 49], under which the three would-be Goldstone bosons {tt"} [50] can be 
formulated as 

U = exp [ivr^r'^/i;] (2.6) 
which transforms, under the SM gauge group S = SU{2)l (g) U{1)y as, 

U ^ U' = g^Ugl , (2.7) 

'^For models where neutrino mass generation is independent of the EWSB sector, see, e.g., Refs. [45, 46]. 
*Mqut stands for the energy scale of the grand unification (GUT) [27]. 
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where gj^ = exp[— i^?£r''/2] and gy = exp[— . Thus the weak gauge boson mass terms, 
My^W^^W^ + ^M'^Z^'Z^^ are given by the dimension-2 nonhnear operator, 

„2 



= ^Tr 



{D^\J){p^V)^ , (2.8) 



where D^U = d^U + z^WI^t^'U - i^UB^r"^ . 

Now consider a pair of generic SM fermions (/, /') which form a left-handed SU{2)l weak 
doublet Fl = (/l, f'jj'^ (with hypercharge Yl) and two right-handed weak singlets fn and /jj (with 
hypercharges Yr and Y^). The electric charges of (/, /') are given by Qf = ^ + Yl = Yr and 
Qfi = —^+Yl = Y^, respectively. We can write down the following gauge-invariant nonlinear 
operator of dimension-3. 



-mfFLu(^l^fR-mf>FLu(^^^^f'n + li.c., 



(2.9) 



which contains the fermion bare mass-terms —mfff — rufiff . 

The SM has no (light) right-handed singlet neutrino z^^, so that the tiny masses for the ac- 
tive neutrinos {vl) indicated by the neutrino oscillations can be naturally of Majorana-type, i.e., 

~nf^^v ^Li^^Li + H-c. , which can be formulated in a gauge-invariant manner, 
2 

= -lei.^Lf^aL^r'$'''e""'e^'3' +H.C. (2.10) 
2 A ■' 

where we choose F^j = Lj (j = 1, 2, 3) as the left-handed lepton doublet, C = ■87^7° is the charge- 
conjugation operator, and ^ = U (o, v/V2j = ■ Thus, the tiny neutrino Majorana masses 
are given by 

= e,,^ (2.11) 

and naturally suppressed by iP/A <^ v in the seesaw mechanism [20] or by the loop-induced fac- 
tor Cij <^ 1 (which also includes small ratios of lepton masses to the TeV scale) in the radiative 
mechanism [21]. If $ is replaced by the linear SM-Higgs-doublet H, Eq. (2.10) recovers the tradi- 
tional dimension-5 operator of Weinberg [19], implying the assumption that the Higgs doublet H 
(responsible for the EWSB) already participates in the neutrino mass generation. However, like 
other SM fermions, the mass generation for neutrinos needs not originate from the same Higgs 
doublet as EWSB [45, 46], so we will not make the ad hoc assumption of <t = H . Instead, we 
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will systematically investigate the most generic neutrino mass operator (2.10) which, in the unitary 

I J, 

gauge, is just the bare Majorana mass term — — Uj^m,yCuj^ + H.c. 

We see that despite the fact that the SM gauge symmetry forbids all bare mass terms for the 
weak gauge bosons and fcrmions, they can nonetheless be formulated in a gauge-invariant way via 
the nonlinear realization and thus are necessarily nonrenormalizable. Hence, the S-matrix elements 
of the high energy scattering involving these massive particles will unavoidably violate unitarity at 
a certain energy scale at which the new physics responsible for the mass generation has to recover 
the proper unitarization. Gcnerically, wc can define the scale A^; for generating a mass nix to be 
the minimal energy above which the bare mass term for rrix has to be replaced by a renormalizable 
interaction (adding at least one new physical degree of freedom to the experimentally observed 
particle spectrum). The unitarity violation scale E* provides an upper bound on such a scale A-r , 
i.e., 

K ^ E\ (2.12) 

and is thus a conservative and model-independent estimate of the mass generation scale Aj, . 

Generally speaking, the scale for EWSB and mass generation of jZ^ is naturally around 
0{v) [cf. eq. (1.1)], but the mass generation scales for the SM fermions are expected to be much 
higher (except, perhaps, for the top quark). To study the unitarity bound on the fermion mass 
we can classify the theories into two types, depending on the nature of EWSB: (i) EWSB is gen- 
erated in a weakly coupled scenario with a fundamental Higgs boson(s); (ii) EWSB is induced by 
dynamical symmetry breaking (with a composite Higgs scalar (s) or no Higgs boson). In either 
case, we can study the scale of the fermion mass generation by examining how the (nonrenormal- 
izable) bare fermion mass term leads to unitarity violation in high energy scattering. For type-(i), 
without assuming the EWSB Higgs boson(s) to also couple to the fermions, we can derive the 
unitarity violation bounds from the scattering // nV^ (which is well-defined above 0(1) TeV 
and below the unitarity violation scale). For type-(ii), due to the dynamical EWSB, the longitu- 
dinal gauge boson Vj^ and its Goldstone boson tt" are composite and cannot be treated as local 
fields at energy scales » 0(1) TeV. Thus, the scattering // nVl no longer provides a reliable, 
model-independent bound on the scales of fermion mass generation. Instead, we may consider the 
scattering // — >■ (//)"" from effective higher dimensional multi-fermion operators. 
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2.2. High Energy Scattering and Power Counting of Energy 

According to the systematic power counting analysis for electroweak gauge theories [51, 52], high 
energy scattering of longitudinally polarized weak gauge bosons exhibits the worst high energy 
behavior and thus is expected to provide the strongest unitarity bounds. This is easy to understand 
intuitively by noting that in the -S-matrix element each external fermion field contributes an energy 
factor oi E2 while an external scalar field gives no additional energy enhancement. However, for 
scattering involving longitudinally polarized weak bosons, each external field gives rise to a 
factor of due to its longitudinal polarization vector 

-] / ^ ^ ^\ ht^ / M \ 

m) = (i^I' = + ""'^^^ ' = ° [-y) ' ^^-^^^ 

where Ma = Mw,z and the term Q{Ma/E) is suppressed in the high energy regime E S> Ma- 
The term k^^/Ma leads to a potentially bad high energy behavior for the scattering amplitude and 
is thus the most sensitive to the unitarity violation. But, the actual energy-dependence of the Vl- 
amplitude is much more involved due to the various intricate energy cancellations. Consider, for 
instance, VlVl — VlVl scattering at tree-level. Naive power counting shows that its amplitude has 
both 0{E'^) and 0{E^) individual terms, but the 0{E'^) terms always cancel out due to the generic 
Yang- Mills gauge structure and the 0(£'^) terms can further cancel down to 0{E^) if a physical 
Higgs boson /i° is present.^ Without such a Higgs scalar h'^, the {VV^, Z°) mass terms may be 
formulated as the nonlinear dimension-2 operator (2.8) which is nonrenormalizable and gives rise 
to non-canceled 0{E'^) terms without respecting to unitarity. This puts an upper bound on the 
EWSB scale [35, 37, 39, 41, 42]. 

The key for carrying out a faithful energy power counting for a scattering amplitude with V^'s is 
to note [51, 52] that the actual i?-dependence is given by in the corresponding amplitude involving 
would-be Goldstonc boson vr^'s. These two scattering amplitudes are quantitatively connected via 
the Electroweak Equivalence Theorem (ET)^*^ [36, 37, 39, 56, 57, 58, 59, 60], 
7 [V^' ,---,V^-; $phys] = C„,od 7 [-iTT-^i , • • • , -iTT"" ; $phys] + B , 

n 

B = Y^ [C"^+i • • • C""?' [v"' , ■ • • , u"" , -Z7r"^+i , • • • , -ztt"" ; $phys] + permutations] (2.14) 

f=i 

= 0{Mw/Ej)— suppressed, 

^It was shown recently [53, 54, 55] that in the compactified (or deconstructed) higher dimensional Higgsless Yang- 
Mills theories, such O(S^) terms are canceled (or suppressed) due to the presence of spin-1 Kaluza- Klein (KK) gauge 
bosons from a geometric Higgs mechanism. 

^"A KK equivalence theorem (KK-ET) [53, 54] has recently been constructed for the geometric Higgs mechanism 
in compactified higher dimensional Yang- Mills theories on general SD-backgrounds (including warped space). This is 
valid independent of whether the zero- mode gauge bosons are massless or not. 
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where v°' = v^V^ and ^phys denotes possible amputated external physical fields [58, 51]^^. The 
B-term on the RHS of the first equation is suppressed by 0(AI\;\/ / Ej) relative to the leading term 
and its precise form [58, 51] is given by the second equation in (2.14). The constant modification 
factor Cmod = C*"^ • • • C*"" = 1 + O(loop) is generated from radiative corrections and can be exactly 
simplified to unity under realistic rcnormalization schemes [57, 58, 59], but will not affect the power 
counting analysis. Since the Goldstone fields vr" have no extra polarization vector like (2.13), there 
is no non-trivial ^'-cancellation in the 7r"-amplitude [the right-hand side of Eq. (2.14)] so that the 
i?-power counting can directly apply.^^ It is shown [51, 52] that such a power counting method 
can be systematically constructed (a la Weinberg [48]) for any given scattering amplitude T and its 
E-powei dependence is, for E 3> Mw^z, fnf , 

Ve = 2 + 2L-^V,Y2-d,-M -e,, (2.15) 
j 

where Vj is the number of vertices of type j which contain dj (^ 0) derivatives and fj (> 0) fermionic 
lines, L is the number of loops, and e„ (^ 0) is the number of external gauge fields v"' = V^v'^ 
appeared in an S-matrix element [cf. (2.13) for definition of v^^]. Prom (2.15), it can be proven [51, 52] 
that, for any n ^ 2, 

De < 2 + 2L , (2.16) 

where the equality holds when the amplitude contains only Goldstone boson vertices. 

In the following sections, we will systematically analyze unitarity for the 2 — > n (n > 2) 
scattering amplitudes ff,U]^i/j^ — nV^ and V^^V^"^ — nV^ . According to the ET discussed 
above, to extract the leading energy behavior, we can analyze the corresponding Goldstone boson 
amplitudes ff,Uj^i/j^ — n7r" and Tr^^Tr^^ _^ n7r" , where the energy power counting can be safely 
carried out [cf. Eq. (2.15)]. At the tree level (L = 0), using the rule (2.15), and the interactions 
(2.9) (which shows that the coupling of two fermions to any number of tt" has no derivatives) and 
(2.8) (which shows that the self-coupling of tt" 's always has two derivatives), we can thus deduce, 

^^This ET relation is deeply rooted in the underlying Higgs mechanism with spontaneous gauge symmetry break- 
ing, even if a fundamental Higgs boson does not exists. As explicitly shown in the second paper of Ref. [57] 
[cf. its Eqs. (88)-(95)] and in Ref. [52], the ET (2.14) is derived from an amputated Slavnov- Taylor identity 
T [Vg^ + iC"i tt"! , • • • , Vg"- + iC"" 7r"" , ^phys] = , which directly reflects the fact that the would-be Goldstone boson 
7r" and the unphysical scalar component of the vector field Vg = egV^ (eg — k'^/Ma, Ma = Mw,z) are "confined" 
together so that no net effect from them can be observed in the physical S-matrix elements a quantitative formulation 
of the general Higgs mechanism at the S-matrix level (independent of whether a Higgs boson actually exists). 

^^In the case that the graviton or its Kaluza-Klein states interact with the longitudinal V^'s via the energy- 
momentum tensor [61], the i?-canccIlation associated with the V^-polarizations again manifests in the corresponding 
amplitude involving the Goldstone boson 7r"'s due to the ET (2.14). 
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0{E') 0{E^) 
(d) (e) 

Figure 1: Typical contributions to // — > nvr'* scattering: (a) the leading "contact" diagram of 
Q{E^)] (b) the sub-leading "t(n)-channel" type diagrams of 0{E^) which can be ignored for the 
unitarity analysis; (c) the "hybrid" contact diagrams (n ^ 3) including Goldstone self-interactions 
from the EWSB sector; (d) the "hybrid" contact diagrams (n ^ 3) including Goldstone interactions 
with the EWSB quanta (illustrated here by the SM Higgs boson as the simplest example of 
EWSB); (e) the "hybrid" contact diagrams due to the assumption of Yukawa interactions between 
the fermions and the SM Higgs boson hP , as the simplest example of fermion mass generation. 

for E » Mw,z,mj , 

7[Vl'VI^ ^ nV£] ~ T[7r"i7r"2 ^ nvr"] 

where the dimensionality of T is always given by D'j = 4— (2 + n) = 2 — n and is independent of 
whether the external fields are bosonic or fermionic. [It will be shown from the exact calculation in 
Sec. 6 that the leading contribution to the second amplitude in Eq. (2.17) is nonvanishing only for 
n(even) = 2,4, 6, • • •.] Eq. (2.17) shows that the E'-power dependence of T is actually independent 
of the number of particles in the final state. Therefore, we note that any conceptual change from 



0(1) 
0(11 



m 



E^ 



E, 



(2.17) 
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including n-body final state reactions (n > 2) has to arise from the n-body phase space integration 
which will be systematically analyzed in the following sections. 

In Fig. 1, we illustrate the relevant contributions to the scattering // m:"'. Fig. 1(a) is 
the leading contribution of 0{E^) from the contact fermion-Goldstone-boson interaction, solely 
reflecting the nonzero bare fermion mass term in the unitary gauge. Fig. 1(b) is the sub- leading 
t/u-channcl type of contributions of 0{E^)^ and can be ignored for deriving the unitarity bound. 
Figs, l(c-e) denote the "hybrid" contact graphs due to the presence of Goldstone self-interactions 
from Eq. (2.8) and the exchange of associated EWSB quanta (such as the SM Higgs). The EWSB 
sector is unitarized at the lower scale of 0(1) TeV, but such contributions in Fig. 1(c) and Fig. l(d-e) 
can be at most of 0{E'^) in general, as we counted above in Eq. (2.17). Unlike the model-indcpcndcnt 
leading "contact" contributions in Fig. 1(a), the additions in Fig. l(c-c) rely on the details of EWSB 
and are thus highly model- dependent. But, this does not affect our general estimate in Sec. 3 which 
only makes use of the above generic power counting results in Eq. (2.17). For a more precise 
analysis, we will first estimate the unitarity bounds by quantitatively computing the universal 
modcl-indcpcndcnt contribution of Fig. 1(a) in Sec. 4.1. Then, in Sec. 4.2, we further examine the 
effect of including the type of contributions in Figs, l(c-e) by specifying an EWSB sector. We will 
show that such model-dependent contributions only affect the unitarity bound by factors of 0(1), 
so that our unitarity limits on the scales of fermion mass generation are robust. 

2.3. General Unitarity Condition and the Customary 2 — > 2 Limits 

For the convenience of later analysis and discussion, we give a unified derivation of the unitarity 
condition on elastic and inelastic scattering including general spins for the particles and a consistent 
treatment of identical particles in the final state, which extends the derivation of Ref. [44]. Then 
we briefly review the customary unitarity limits for various elastic and inelastic 2 — 2 scattering 
processes. 

2.3.1. General 2 — ^ n Unitarity Condition 

The unitarity bound originates from the unitarity condition of the S-matrix, S^S = 1, which, with 
the definition § = 1 -|- z T , can be expressed as, 7^7 =2 SJm 7 . Taking the matrix element of both 
sides of this relation between identical 2-body states and inserting a complete set of intermediate 
states into its left-hand side, one arrives at 




(2.18) 
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where on the left-hand side Jpg denotes the n-body phase space integration (including the sym- 
metric factors from the possible identical particles in the final state) and, in the second term, the 
summation is over all inelastic channels (n ^ 2). The right-hand-side is evaluated in the forward 
direction. For the 2 — >■ 2 elastic channel, we make the partial wave expansion for the scattering 



amplitude 



Tei[2^2] = 167re^(^-^>^(2j + l)<,^(cose) 



af 



1 .1 ' _ (2.19) 



el _ ^ i(v'~u)ip 



e'^"" '''^'f / dcosO <,^(cos e) Tei[2 ^ 2] , 



^ 327r 

where 6 and are the scattering angles, and [v, u') = {vi — V2^ '^3 — '^4.) with (ui, U2) and (^3, V4) 

the helicity indices of the incoming and outgoing particles, respectively. The d-functions satisfy the 

1 2 (5 •/ 

dx dl,,^{x) dl,,^{x) = '''' , and for the special case of v = u' = , we 
-1 2_7 -|- 1 

have dQQ{x) = Pj{x) , with Pj{x) the Legendre polynomials. So, for the incoming and outgoing 
particles of the same helicities, ui = 1/2 and 1^3 = 1^4, (2.19) reduces to 



Tei[2^2] = 167^^(2j + l)P,•(cos^)a^^ 

T .1 ' (2.20) 
af = — dcos0Pj(cos0)Tei[2^2]. 

OZTT J_i 

With (2.19), it is straightforward to compute 

/ |Tei[2^2]|2 = ^Y.(^j + l)\af\\ (2.21) 

JPS2 J 

where the symmetry factor equals 1! (2!) if the 2-body final state consists of non-identical (iden- 
tical) particles-*^^ . In (2.21) and the rest of this paper, we consider the mass of any initial/final state 
particle to be much smaller than the cm. energy, rOj <C s , which is justified for the scattering 
processes to be analyzed in Sec. 3-6.^'^ Using Eqs. (2.19) and (2.21), the unitarity condition (2.18) 
becomes 

+ 1)1 [I - [^cafY - {<^maf - = -L ^ /" 17,^^,[2 n]f ^ , (2.22) 



Because the right-hand side of (2.22) is non-negative, we find a consistent solution, for each j , 



^taf) +{Qmaf-^) ^ |, (2.23) 



^^Throughout our derivation the factors due to the proper counting of final state identical particles will solely 
originate from the phase space integration of the n-body final state in the total cross section. 

'^"'For completeness, in AppendixB we give the precise formula of /pg^ |Tei[2 — » 2]f and the refined unitarity condi- 
tions by including the possible masses of initial/final state. 



14 



which results in 



el 



of 



Qe_ 

2 ' 

Be- 



(2.24) 



When the 2-body final state consists of non-identical particles [q^ = 1), these conditions reduce to 
the familiar form, |3fiea|'| ^ ^ and \af\ ^ 1 . With (2.24) and defining the elastic cross section 

(Tel = ^ af , we deduce a bound on cTei in its j-th partial wave, from Eq. (2.21), 
j 

167r(2i + 1)^, 



<[2^2] ^ 



(2.25) 



which, for = Ij agrees with Eq. (38.47) of Ref. [62]. Finally, by assuming the dominance of 
j = contribution to the elastic channel, an upper bound on the inelastic cross section crinei[2 n] 
(n > 2) can be derived from Eq. (2.22), 



CinelP 



(2.26) 



Here the identical particle factor Qg is determined by the final state of the elastic channel (which 
has the same initial state as the inelastic channel) as indicated by the subscript of Qe ■ For a given 
initial state of the inelastic channel, as long as we can find at least one elastic channel (with the 
same initial state) in which the final state particles are non-identical^^ , we can then set Qe = ^ in 
the above inelastic bound (2.26) to get the optimal constraint. This mild requirement is satisfied 
for all the inelastic channels of //, //' nV^ . Before concluding this subsection, let us derive a 
partial wave unitarity condition for the 2 — 2 inelastic channel [parallel to (2.24)]. Using the same 
partial wave expansion formula (2.20) for the inelastic amplitude, we can compute 

327r 



/ 



|Tinel[2^2]r 



E(2i + i)l«: 



inel|2 
'3 I ' 



(2.27) 



where q-^ = l! (2!) for an inelastic final state being non-identical (identical) particles. Substituting 
this back into the RHS of (2.22), we deduce 



> 



V(2j + l)-|af<=i'2 



3 



1 



(2.28) 



so that, for each j , 



inel I 



< 



(2.29) 



^^For some types of initial states in a given model, the final state with identical particles may be automatically 
forbidden due to the particular charge of this initial state or the absence of certain couplings. 
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Interestingly, the optimal bound for a given 2 — >^ 2 inelastic channel is realized when we can set 

Qe = I for a corresponding elastic channel (which shares the same initial state), i.e., |a™*^^| < . 
We see that this or (2.29) takes a different form from the elastic partial wave condition (2.24) when 
final state consists of identical particles or when aj . 



2.3.2. Customary Limits from 2 — ^2 Scattering 

2.3.2A. 2^2 Unitarity Bound for the EWSB Scale 

The unitarity violation from the quasi-elastic scattering V^^V^ — V^^V^* was originally 
studied by Dicus-Mathur/Lee-Quigg-Thacker (DM/LQT) [35, 37] for constraining the Higgs boson 
mass in the SM. According to the ET, its leading high energy behavior is represented by the 
corresponding Goldstone boson scattering Tr^^Tr^^ — ^ tt^^tt""^ and can be directly computed from 
the Higgsless nonlinear Lagrangian (2.8) as 

TfTrOTrO ^ 7r\^] = , T[7r07rO ^ Tr+Tr"] = ^ , T[^+^- ^ ttV"] = A. , 

t s u (2-30) 

where t ~ —^s (1 — cos^) and u 2± —^s (1 + cosO) . This is usually called the low energy theo- 
rem [63] and may be reexpressed in terms of isospin amplitudes 7[I] [64] , 

T[0] = 3A(s,t,u) +A(t,s,u) +A(u,t,s) = ^ , 

t — u 

T[l] = A{t,s,u)-A{u,t,s) = — (2.31) 



T[2] = A{t,s,u) +A{u,t,s) 



t + u 



g 

where Ais, t, u) = , and the relevant identical particle factors arc not included. Considering 



the isospin-singlet channel and including the identical particle factor = 21 in the first inequality 
of Eq. (2.24), we see that the unitarity condition |3?eao| ^ ^ leads to an optimal limit, ^/s ^ 
with = \/Sttv ~ 1.2TeV. Another derivation is to define the normalized isospin-singlet state 
(including identical particle factors), |0) = — ^ [2 |7r"'"7r~) -|- Itt^tt'^)] , and compute the scattering 
amplitude in the singlet-channel, 

T[0] = \ [4T[+-, +-] + 4T[+-, 00] + T[00, 00]] = ■ (2.32) 

Thus, from (2.32) the usual unitarity condition of the s-wave amplitude ( |3fteao| ^ 1/2) imposes 
an optimal bound for the isospin singlet-channel, 

= V^v ~ 1.2 TeV. (2.33) 
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which we mentioned ear Her in Eq. (1.1). 



2.3.2B. 2 — > 2 Unitarity Bound on the Scale of Fermion Mass Generation 

Appelquist and Chanowitz (AC) [43] considered the scattering f±f± — V^V^ [65] with fermion 
bare mass terms^^, and derived the asymptotic amphtude, 

7\!±]±^WlWl\ = ±^V~s, (2.34) 

where the subscripts for f±f± indicate the fermion hehcity zt— . For the color-singlet initial state, 

this gives the s-wave unitarity limit s/s ^ — , as we mentioned earlier in Eq. (1.2). 

VNcirif 

Defining the spin-singlet combination for the initial state fermions, — -= [|/_|-/+) — |/_/_)] , we may 

v2 

make this bound slightly tighter, 

E*f ~ ^= . (2.35) 

One can further define a normalized isospin-singlet state for the out-going gauge bosons with the 

combination [41] ^ [2|W2"^r) + I-^l-^l)] > and thus slightly reduce the bound (2.35) by an 
v6 

additional factor of \ I ~ 1/1.22 , i.e.,^'^ ~ j!Z . 

V 3 ^ ' ' f ^/SWcmf 

2.3.2C. 2—^2 Unitarity Bound on the Scale of Majorana Neutrino Mass Generation 

The Majorana neutrino scattering I'j^^i^j^^ ^l^^l^ '^^^ considered recently by Maltoni- 
Niczyporuk-Willenbrock (MNW) [44, 66] who derived an upper bound on the scale of the neutrino 
mass generation by using the Weinberg dimension-5 operator [19] with the usual SM Higgs doublet. 
For the inelastic scattering — or -^{hh), one obtains, from 

the condition (2.29) (^>e = 2),^^ 

E: . i^. (2.36) 
m,. 



^^Thcir equivalent gauge- invariant nonlinear formulation is given in Eq. (2.9). 

^^For clarity of the following analyses with a multiple V£ (tt") final state, we will ignore such minor improvements 
and thus will use (2.35) for comparison. 

^*In this paper we will not add extra identical particle normalization factors for the initial state as they are irrelevant 
to the calculation of the cross section and unitarity bound (cf. Sec. 2.3.1). For convenience, we may do so in the 
amplitude calculation only for final state identical particles as far as the total cross section and unitarity bound are 
concerned (which also means that in the condition (2.29) we will accordingly remove the final state identical particle 
factor and still retain Qe), but we keep in mind that this is inappropriate when computing the differential cross 
section. 
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which can be improved to ~ 



by using a mixed final state \\ZlZl + hh) [44]. 



In our analysis, the bare neutrino Majorana mass term arises from the nonlinear operator 
(2.10) without assuming the participation of the SM Higgs boson. With this, the scattering 



processes 



V2 



'1W^).^\ZlZl) and 



Itt+tt ), ;^|7r°7r'^) can be computed. The helicity amplitudes are 



TT TT 



~ ±^r-^/s , 



7 



^l±vl± 



V2 



ZlZl 



~ -7 



V2 



-^0 
TT TT 



Thus, for a spin-singlet initial state, the scattering 
to a unitarity violation limit 



\^L+^L+) 



E* ~ 



2V2 



TTV 



(2.37) 

. -^\ZlZl) leads 
(2.38) 



The reason that this limit is stronger than (2.36) by a factor of 2 is the absence of the SM Higgs 

boson in our nonlinear operator (2.10). This is due to the fact that the s-channel SM-Higgs exchange 

arising from the usual dimension-5 Weinberg operator will partially cancel with the gauge amplitude 

in the unitary gauge and push the unitarity violation to a slightly higher scale. We also note that 

with the isospin-singlet combination for the final state gauge bosons, [2\W^W£) + IZ^Z^^)] , 

V6 



the bound (2.38) may be further improved by factor of 
mass rrij, ~ 0.05 eV, the limit (2.38) [or (2.36)] gives, E^ 
unification (GUT) scale. 



V3 



1/1.15 . With the typical neutrino 



10 GeV, which is right at the grand 
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3. Challenge from 2 ^ n Inelastic Scattering 

3.1. Puzzle: Energy Power Counting vs. Kinematic Condition 

According to the power counting result Eq. (2.17), the E-power dependence of the scattering am- 
phtudes for fj.v^v^ nVf; and V^^V^^ nVf; (n ^ 2) is independent of the number of V^'s 
in the final state. This means that any conceptual change for n > 2 has to come from the n-body 
phase space. Indeed, it is straightforward to count the £'-power dependence from the n-body phase 
space integration. 



— f 

3in Jp 



^ 1 r ^-Uv^-X ^^ ^4-3n^(4) 

PS„ gSinJ 2E,...2En^ ' / (3.39) 




^ ^ £;2(n-3)^ 



where 0in = 4 [(pi • — (mim2)^] ^ ~ 2s is the flux factor of the initial state, q is the symmetry 
factor from possible identical particles in the final state, and we define E = -s/s . Eq. (3.39) shows 
that the phase space contributes a nontrivial energy power factor to the cross section which grows 
with n. 

Combining the powers of energy for the scattering amplitude in Eq. (2.17) and for the n-body 
phase space integration in Eq. (3.39), we arrive at the following naive power counting estimate, 
including all dimensionful parameters in the cross section, 

a[2^n] cx - (^^j (^) , (n ^ 2) , (3.40) 

where 6 = 1 for ff,u^v^ nV^ (as considered in [44]) and 5 = 2 for V^^V^^ — > nV^ . For 
simplicity, we have used to denote the mass of cither a Dirac fermion / or a Majorana neutrino 
Vj^. With the general condition (2.26) the naive estimate (3.40) results in the unitarity limit. 



E* V 



/ \ 2(2-5)' 



2(n-2+d) 

— > V , (for n — s- large) , (3-41) 



where Cq is a dimcnsionlcss constant undetermined by the power counting. From (3.41), wc sec that 
for 6 = 1, it gives the bound (1.3) (which is the conclusion of Ref. [44]), while for 6 = 2, it results 
in the bound (1.5). It is clear that if the above naive power counting estimate could really hold, 
then one is led to a striking conclusion that all unitarity bounds from multiple V^-production would 
approach the same scale v ~ 246 GeV for arbitrarily large n, and thus such scattering processes 
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would reveal no new scale (other than v) for the mass generation. As we have observed, even for 
the EWSB scale, this would imply a significant reduction of the customary 2^2 limit (1.1) [or 
(2.33)]. 

The puzzle is sharpened when we consider the general kinematic condition for all multiple Vl- 
production processes, 

^Ts > nMw(z) ^ 3 ^ , E* > V-, (3.42) 

which must hold for all these processes as we advertized in the Introduction. When n becomes 
large, this condition requires the unitarity violation scale = ^* to grow at least linearly with 
n , rather than arbitrarily approaching the scale v . This unavoidably contradicts to the conclusion 
drawn from the naive power counting estimate (3.41). It is this observation that leads us to believe 
that the condition (3.41) cannot really hold and a deeper resolution must exist. 

3.2. Resolution: Increasing Power of Energy vs. Phase Space Suppression 

At this point it appears clear that in order to be consistent with the kinematic condition (3.42) there 
must be additional n-dependent dimensionless factors in the exact n-body phase space integration 
that sufficiently suppress the E-powei enhancement in Eq. (3.39). 

When the possible angular dependence is ignored in the squared amplitude, |T|^, the cross 
section factorizes into the squared amplitude times the phase space, so the n-body phase space 
integration can be done exactly (cf. Appendix A) and is given by 



1 f 1 fd^ki---d^kn 



Sin Jps„ Q3 

„n— 3 




(3.43) 



24{«-i)7r2"-3(n- l)!(n-2)!^ ' 

For the processes //, I'li']^ nVl and V^^V^"^ —>■ nV^ , using Eq. (2.17) we can generically write 
the form of the squared amplitudes as, 

= c,{9j) (2N,f-' , (3.44) 

where S is defined below Eq. (3.40) and CQ{9j) is a process-dependent, dimensionless coefficient with 
possible dependence on the scattering angles (Oj)^^. Here, the factor 2Nc comes from the color- 
singlet channel of the initial state fermion pair with a proper helicity combination [cf. Eq. (4.58)], 

^^It turns out that for the scattering //, UlI^l riV^ , the universal leading contributions are given by the "contact" 
Feynman diagrams represented in Fig. 1(a) where the Cq has no angular dependence. 
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and Nc = 3(1) for quarks (leptons or neutrinos). Thus, we can have a realistic estimate for the 
scattering cross section, 



24("-i)7r2"-3(n- l)!(ra-2)! s 

where Cq = Cq/ g and Cq is a constant factor originating from the coefficient Cq in the squared 
amphtude (3.44) and includes the possible modification to the n-body phase space integration 
(3.43) (if Cq has nontrivial angular dependence). Substituting (3.45) into the general condition 
(2.26), we thus arrive at the following improved estimate of the unitarity bound. 



E* 



Co 24"-2;r2(n-i) (n - 1)! (n - 2)! {2Ncf-'^ 



. I V X 2(2-5)- 
5-2 I V 



ruf 



1 



2{n-2+«) 

(3.46) 



where the process-dependent constant Cq will be evaluated by an exact calculation (cf . Sec. 4-5) . 
For the current estimate it is reasonable to assume Co = 0(1) or, in a more prescise form. 



which becomes increasingly more accurate as n gets larger. In fact, Eq. (3.46) may be regarded as an 
exact bound (with the constant Cq to be derived from the precise calculation), and thus the formula 
(3.47) may actually represent all our approximations made in the current estimate of the unitarity 
bound: (i). ignore the dimensionless coefficient Co(%) in the amplitude-square (3.44) including 
its possible angular dependence [that may lead to an 0(1) correction to the n-body phase space 
integration (3.43)]; (ii). ignore the symmetry factor l/g arising from possible identical particles in 
a given final state. 

The crucial observation from Eq. (3.46) is to note that the n-dependent dimensionless factor 

1 

[(n — 1)! (n — 2)!] 2(n-2+d) jQgg approach unity as n becomes large, rather, it increases almost 
linearly with n. This can be understood from Stirling's formula^*^, n! ~ n"e~"\/27m, from which 
we can estimate, for large n, 

1 1 77 72 

[(n- l)!(n-2)!]2("-2+«) ~ (n!)n ~ - > -, (forn»l), (3.48) 

e 3 

where the constant e = 2.7182818- • •. It is this feature that makes our bound (3.46) fully consis- 
tent with the kinematic condition (3.42)! For the processes ff,Uj^Uj^ nV£ {5 = 1) with light 
quarks (leptons) or neutrinos, i.e., f^/m^ ^ 1, we observe that there is a competition in the uni- 
tarity bound (3.46) between the quickly decreasing factor [y / m "--'^ and the almost linearly 
^"The relative error of this formula is less than 4% for n > 2 . 
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increasing factor [(n — 1)! (n — 2)!] 2{"-2+5) when n becomes large. As a consequence we expect 

the strongest unitarity bound to occur at a value n = Ug > 2 and, after n exceeds Ug , the factor 

1 

[(n — 1)! {n — 2)!] 2{n-2+5) become dominant and eventually force the bound (3.46) to grow 

almost linearly with n . Only for the heaviest SM fermion, the top quark (where v'^ /ml ~ 2 ^ 1), 
can the value of Us be equal (or very close) to 2. On the other hand, for the process V^^V^'^ nV^ 
((5 = 2), when n increases, the only possible decreasing factor on the right-hand side of Eq. (3.46) 
is (^^/Co / (27r)^ " , provided y/Co > 2vr . But our current estimate under the approximation (3.47) 
implies the setting Co ~ 1 . So it seems that we will have = 2 for the EWSB scale which 
corrsponds to the customary 2^2 unitarity hmit discussed in Sec. 2.3.2 [cf. Eq. (2.33)]. Since the 
current estimate does not precisely fix the value of Cq, a reliable determination of Cq and thus 

should be given by a quantitative calculation, which will indeed support Ug = 2 (cf. Sec. 6). 

1 j_ 

Considering the region n S> 1 and thus 0(1) 2("-2+'5) ~ 0(1)2" ~ 1 , we can further derive a 
complete asymptotic formula for Eq. (3.46), 



E* 



47r 



m, 



Attv 



-n , (for n S> 1 ) 

(2 - 5)2 



n + (2 - (5) In + 



m 



n 



(3.49) 



Prom this, we deduce that for each asymptotic value of n, the leading difference between the unitarity 
bounds of two types of fermions with masses m^^ and m^^ is given by^-*^ 



El 



Attv 



In 



m. 







vrii 



n 



(3.50) 



which is independent of the value of n. Similarly, the difference between the unitarity limit E*^ for 
the fermion / and the limit E^ for the scale of EWSB behaves as 



Airv 



In 



+ 



m 



n 



(3.51) 



and is again independent of n in its asymptotic region. 

Finally, in Fig. 2, we numerically plot out the bound (3.46) under the approximation (3.47), for 
a number of typical processes including V^^V^"^, ti, bb, t~t^, e~e+, I'l^i'i —>■ nVl , where we have 
^^For the analysis in the current section, / denotes either a quark, a charged lepton, or a Majorana neutrino. 
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input the following central values for the SM fermion masses [62] 

{mt, rub, rUc, rus, ruu, ma) = (178.0, 4.85, 1.65, 0.105, 0.003, 0.006) GeV, 

(3.52) 

{mr,m^,me) = (1.777, 0.1057, 0.000511) GeV. 

For leptons and heavy quarks (c, b, t) the pole-masses are used, while for light quarks (n, d, s) the 
current-quark masses under MS [62] are chosen as the light quark pole-masses are hard to directly 
extract. Our 2 — n analysis shows that for light fermions the minimum n = Ug is sizable and the 
unitarity limit is insensitive to the fermion masses, so a small difference between the pole-mass and 
MS mass for light fermions does not cause a visible effect in Fig. 2. The absolute scale for neutrino 
masses will be further determined by various experiments. As mentioned earlier, the recent WMAP 
data [14], in conjunction with the 2dF Galaxy Redshift Survey [15] put a 95%C.L. upper limit 

^TTT-i/j ^ 1.01 eV [16]. From the WMAP, SDSS and Lyman-a Forest data, a stronger bound [17], 

j 

mnj ^ 0.65 eV (95%C.L.), may be derived. Neutrino oscillations measure the difference of 
squared mass-eigenvalues and the atmospheric data require the larger mass gap bounded by [4, 10], 

0.04 eV ^ Aft^ ^ 0.06 eV, (99%C.L.), (3.53) 

where A^tm = !"^^i2 ~ '^vj • This requires the neutrino mass scale to be bounded from below, 

i.e., rrii, > A^^^^^ ~ 0.05 eV, and when the neutrino mass spectrum exhibits a hierarchical structure, 
1 

the quantity A^^.j^ actually sets the mass scale for one or two of the active neutrinos [10, 69]. 
The upcoming laboratory experiments on neutrinoless double /3-decay [70, 71] are indispensable to 
further pin down the neutrino mass scale. Without losing generality, wc choose a sample neutrino 
mass value rriv — 0.05 eV for the analyses in Fig. 2 (and in the Table 1 below). 

Fig. 2 shows that the scattering processes //, Vifi nV^ indeed put new unitarity limits on 
the scale of fermion mass generation, independent of the scale v, and in agreement with the original 
motivation of Appelquist and Chanowitz [43]. However, strikingly, our analysis further reveals that 
for the light fermions including neutrinos, leptons and quarks, the strongest bounds generally occur 
at a value n = Ug > 2 . The current improved estimates of (n^, E^^^) are summarized in Table 1. 
Note that these bounds easily satisfy the kinematic constraint Eq. (1.4). 

Another feature of Fig. 2 is that the curves for the quarks and charged leptons become almost 

parallel in the asymptotic region n > 15 , and the upper curve (for neutrinos) also becomes parallel 

when n > 35. This fact can be understood from the asymptotic formula, Eqs. (3.50) and (3.51). 

^^Rcccntly both DO and CDF announced their updated analyses of the Run-I data for top quark mass measure- 
ments [67]. The resulting world average for the top quark pole-mass is mt = 178.0 ± 4.3 GeV [68]. 
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Figure 2: Realistic estimate of the unitarity bound E*, based upon Eq. (3.46) with the approximation 
(3.47), for the scattering processes V^^V^'^, tt, hh, t~t~^, e~e~^, Vj^V]^ — > nV^ (curves from bottom 
to top), as a function of n (^ 2). Only those points corresponding to integer n have physical 
meaning, and on the lowest curve n is restricted to the even integers. 

Fig. 2 further shows that for the scattering tt nV£ the unitarity bounds are very close for 
n = 2,3, 4, and it is easy to check that the minimum of the bound can shift to = 4 if we 
just average over the spin and color for the ti pair instead of composing a spin-0 and color-singlet 
initial state. The precise calculations in Sections 4, 5, and 6 below will quantitatively improve 
the above estimated numbers but, as we will see, all the qualitative features in Fig. 2 and Table 1 
remain. Finally, it is very useful to compare these estimated bounds (Table 1) with the customary 
2 — > 2 limits (2.33), (2.35) and (2.36) [cf. Sec. 2.3.2] whose predictions are summarized in Table 2. 
Comparing Table 1 and 2, we see that for the scales of EWSB and top mass generation, our estimates 
in Table 1 agree with the customary limits in Table 2; however, for all the light fermions and for 
the Major ana neutrinos, our bounds are substantially stronger than the AC [43] and MNW [44, 66] 
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limits by many orders of magnitude. 



Table 1: The estimated lowest unitarity bound E^^^ [derived from the Eqs. (3.46) and (3.47)] for 
the scattering process ^^^^2 ~^ "-^l '^^^^^ the corresponding number of final state particles n = Ug , 
where ^12 = Vl, f, or their anti-particles. 



^1^2 




ti 


bb 


cc 


ss 


dd 


uu 




11-11+ 


e 6+ 




Us 


2 


2 


6 


7 


10 


12 


13 


7 


10 


15 


32 


^min(TeV) 


1.2 


3.5 


14 


18 


26 


35 


37 


19 


28 


44 


95 



Table 2: Predictions of the customary unitarity limits for the scattering process ^^^^2 ~^ ^l^^l^ > 
based upon Eqs. (2.33), (2.35) and (2.38). 
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4. Scales of Mass Generation for Quarks and Leptons 



In this section we will attempt to improve on our new limits on the scales of fermion mass generation 
by computing the matrix elements more precisely. A parallel treatment of Majorana neutrinos will 
be given in Sec. 5. In Sec. 4.1, we first systematically derive the precise unitarity limits for various 
// — > nvr" {nV£) scatterings with quarks/leptons in the initial state, which result in optimal upper 
bounds on the scales of the fermion mass generation. Then, we compare these quantitative limits 
with our estimates given in the previous section (Table 1 and Fig. 2). In Sec. 4.2 wc include EWSB 
effects among the vr" [cf. Fig. l(c-d)]. We also comment on fermion mass generation with dynamical 
EWSB, where the V£ and its would-be Goldstone boson tt" are composite fields. 

4.1. Uniteirity Bound on the Mass Generation Scale from // — > riTr" 



int 



-T- 



n=l 



v"'n\ 



Prom the nonlinear fermion mass term (2.9), we deduce the following fermion-Goldstone interaction 
Lagrangian, 

|7f|" {rrifff + mf,'f'f'^ , {n = even) , 

+ V^TT-f (mfPL - mj,PR^ f + H.c] , (n = odd) , 

where |7f| = [27r"'"7r~ + tt^tt'^] ^ and Pl{R) = ^(1 T 7^)- This leads to three types of 2 — > n 
processes, containing the leading contact diagrams shown in Fig. 1(a), 

(a) . //, f'7 - K)' {n-y i^T'' , (n = even) , 

(b) . //, f'T - K)' (^-)' i^T'"' ^ iji = odd) , (4.55) 

(c) . ff - {n'-y^' {n-y {^r~"'~' , = odd) , 

which may be generically denoted as /"/^ (7?+)^^ (vr")^ (vr^)" ^ ^, where /",/'' G (/, /') as 
appropriate. The helicity amplitudes T f^f^_^;n,k,£ are given by 



n 



r — 1 nifs/s 

even: T[/±/±; n,i,i] = ^i^Ci^^ 



n = odd: T[/±/^;n,£,£] = -i"e2^^^, 



J 



f±f±;n,e,e 



/±/l;n,£ + l,^ =-i''^^[Z! a'[/4/±;n,£,^ + l] = -i-e3(7 



nif, 



rrif J ti" ' 



(4.56) 
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where the coefHcients (Ci, 62, C3) are defined as 



60 = 



e. = 



1 

n! 2 
1 



CLi2^ i£\f{n-2£)\, 



(4.57) 



^-CLi2U\ {£ + iy. (n-2£- 1)!, 



n! 



and ci = 



Now we define a color-singlet initial state with a proper helicity combination, 



|in) 



1 ^= n 



/3=1 



ja/3 



(4.58) 



/3=1 



where on the right-hand side Nc = 3(1) for quarks (leptons), and for the first definition the sign 
— (-I-) between the two helicity states corresponds to the n = even (odd) case for an initial state 
|in) of spin-O(l). With these we compute the cross section for |in) — {n'^)'^ {t^~)^ ij^^)^ ^ ^ as 



o"[in; n, A;, l\ = 



Coj = — 



1 fVsV^^'-'^f'^f^' 



24("-i)7r2"-3(n- l)!(n- 2)! s \ v 

1 

Q 



p2 

^3 



(4.59) 



k\l\[n-k-l)\ ' 

where / G (/, /'), and x = 1 (0) for the first (second) definition of the initial state in Eq. (4.58). 
For the three processes in Eq. (4.55), the maximal values of Cqj can be derived as 

2" (§!)' 



Cmax 
01 — 



(n!) 



2 ' 



Gmax 
n 



^02 



^max _ ^ V 2 /• V 2 /• 

"7 ' *-^03 — 



(n!)^ 



(4.60) 



Tt Tl -f- 1 

corresponding to k = i = —, k = £ = 0, and k = £ + 1 = — ^ — , respectively. Thus, for each given 
type of fermion, using the general unitarity condition (2.26), we arrive at the bound -^/i ^ E*^ with 

(4.61) 



= V 



V 



jmiax given by 



/^max 



max 
■j 



2(n-l) 



24(n-l)^2n-3(„ 2)! ' 



(4.62) 
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or, using (4.60), 



<T)max 

mm; 



mmax 









23n- 


-4^2n-3 (n!)2(n-l)!(n 


-2)! 




1 




24(n 


-i)7r2«-3n! (n- l)!(n- 


-2)! ' 








23n- 


-47r2"-3 (n!)2(n- l)!(n 


-2)! 



(n = even) , 
(n = odd) , 

(n = odd) . 



(4.63) 



where 0^™!^ have x = 1, 1, and correspond to the processes (a), (b) and (c) in Eq. (4.55), 
respectively. 

Table 3: The precise minimum unitarity bound i^j^j^ [derived from Eq. (4.61)] for the scattering 
process ,^^^2 ~^ ("-^l) corresponding number of final state particles n = Ug , where 

Ci,2 = foif. 



^1^2 


tt 


hh 


cc 


ss dd 


uu 


T T'^ H n'^ e 


Us 


2 


4 


6 


8 10 


10 


6 8 12 


^min(TeV) 


3.49 


23.4 


30.8 


52.1 77.4 


83.6 


33.9 56.3 107 



Finally, using the bound (4.61), we derive the lowest bound for each scattering and the cor- 
responding number of final state particles n = Ug , as shown in Table 3. This is to be compared 
with our estimated bounds in Table 1 and the customary limits in Table 2. We see that the precise 
bounds agree well with our estimate for the case of a heavy top quark, but significantly improve 
the estimates in Table 1 for all other light fermions; the limits in Table 3 become about a factor of 
1.5 — 2 higher. The shape of each bound as a function of n is shown in Fig. 3 for all SM quarks 
and Icptons, where Ug corresponds to the minimum of each curve. The values of Ug are visibly 
smaller than those in Table 1 and Fig. 2. Once n > Ug , the curves also grow much faster than 
Fig. 2, but still exhibit very similar slopes for n > 15 . This is because for the large n the power 
factor (v/nif)"^ is essentually unity and the dimensionless phase-space factor [cf. Eqs. (3.48) and 
(4.62)] becomes dominant. Using the Stirling formula, we find that the three functions 0^™!^ in 
Eq. (4.62) exhibit the same asymptotic behavior 

2n 

(forn»l), (4.64) 



rnmax 
■^1,2,3 



3 
62 
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so that the final unitarity limit (4.61) behaves as 



' asym 
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Att 3 












62 



( for n 3> 1 



Attv 

3 
62 



n2 + n2 In 



1 _l , 2 

+ —n 2 In 



(4.65) 



(4.66) 



where we approximate 0(l)2("--i) 1 for n S> 1 . We see that unlike Eq. (3.49), the leading term 
in the above bound has (rather than n^) power dependence. This is why the curves in Fig. 3 
grow much faster than those in Fig. 2 in the asymptotic region of n . Furthermore, we note that 
the second term in the expansion (4.66) depends on both and m-, contrary to (3.49). But, the 
third term in the expansion of (4.66) is only suppressed by n~2 (not n~^) and thus can be sizable 
due to the squared logarithm. For instance, comparing the limits for 6~6+ and tt scatterings and 
using Eq. (4.65), we obtain, numerically, 

61.1 TeV, (n = 15), 
63.4 TeV, (n = 30) , 



'asym 



(4.67) 



which explains why in Fig. 3 the curves are almost parallel for the displayed asymptotic region 
15 < n ^ 30. Finally, the charm quark has a slightly stronger bound than the tau lepton despite 
nic < rrij- (cf. Fig. 3 and Table 3) because the color factor Nc = 3 in (4.61) helps to lower the 
bound for charm. 
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10 20 30 

n 

Figure 3: Precise unitarity bound E*j based upon Eq. (4.61), for the scattering processes 
ti, bb, cc, T~T^, ss, iJ^'fJ-^, dd, uu, e~e^ — > mr^ ('^^l) (curves from bottom to top), as a func- 
tion of n ( ^ 2 ) . For the points with odd n values, the best bounds are actually given by 
C-I-S+, T^i^-, C-s^, l^-i^-, U-d-, e^u- —>■ mr"- (curves from bottom to top), re- 

spectively. Only integer values of n have physical meaning. The curves marked by circles represent 
bounds for leptons and those marked by triangles are for quarks; also, for clarity, the up-type, 
down- type quarks and the leptons are plotted in three different colors. 



4.2. Model-Dependent Effects from EWSB Sector 

As shown in Fig. l(c,d), for n ^ 3 the tree- level scattering process // mr"' can have additional 
contributions involving the EWSB sector, i.e., the Goldstone boson self-interaction vertices and 
the exchange of the new quanta that unitarize the EWSB sector (such as the SM Higgs boson). 
These contributions depend on the details of the EWSB mechanism and are thus highly model- 
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dependent. In the following, we will take the simplest EWSB sector of the SM as an explicit 

example and examine how our universal estimates of the unitarity bounds [based upon the pure 

contact interactions among fermions and Goldstone bosons in Fig. 1(a)] are affected by the EWSB 

contributions via Fig. l(c,d). The diagrams in Fig. 1(a) and (c,d) are of the same 0{E^) by power 

counting, so that a contribution of type-(c,d) may affect the amplitude by a factor of 0(2 — 3) 

and thus the total cross section by a factor of [0(2 — 3)]^ ^ 0(10) . Thus we expect the unitarity 

1 

bound of Sec. 4.1 to be changed by at most a factor of [0(2 — 3)] ("^i) (n ^ 3) since the n-body 

phase space integration contributes an energy factor of s"^^ (^^f^ Sec. 3.1-3.2). Numerically the 
1 

factor [0(2 — 3)] ("^i) is at most 1.7 (for n = 3), and becomes close to unity as n increases above 3. 
This means that, fortunately, the effect of the EWSB contribution is small for light fermions since 
their strongest unitarity bounds occur at significantly large n values (cf. Fig. 2-3). The effect for 
top-quark case is interesting since the analysis based on contact contribution alone shows that the 
unitarity bound is strongest for = 2 and only becomes slightly weaker for Ug = 3. Hence, adding 
the EWSB contribution associated with tt nvr'^ (n ^ 3) processes could push the location of the 
minimum unitarity limit up to Ug = 3. As explicitly shown below, this is indeed the case, but the 
improvement in the unitarity bound for all fermions (including the top-quark) is always less than a 
factor of 2, as expected. 

For convenience we first consider the SM with its EWSB sector (Higgs sector) nonlinearly real- 
ized. The Higgs-Goldstone fields of the SM can then be formulated as 2 x 2 matrix, 

^ = {v + h?)U , U = exp[i7r''r"/t;] , (4.68) 

where the physical Higgs field hP transforms as a singlet under SU{2)l <8) U{1)y- Correspondingly, 
the EWSB Lagrangian of the SM takes the form, 

^EWSB = - Vih) = \Tr:[{D^^){D^^)^] - \[{v + h^f - v^]\ (4.69) 

where D^^ = {d^h^)U + {v + h^)D^U and D^U is defined below Eq. (2.8). The kinematic term 
can be simplified as 

= ^(aM^°)(5'^/i°) + ^(^^' + 2i;/iO + /t°')Tr[(L»'^i7)(Z?^C/)t]. (4.70) 
The part —Tr[{D^'U){D^U)^] is the same as Eq. (2.8) and contains the following pure Goldstone 
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Lagrangian, 

£gb = ^-^Tr[{d^U){d,Uy] 

(If • di^itf -{If -If) (d'^it ■ d^if)] , 

rt(cvcii)=4 

(4.71) 

where 

It - It = 27r+7r-+7r°7r° , 

( vf • d^ltf -{If -If) {d^lt ■ d^lf) (4.72) 
= (7r+9^7r-+7r-5^7r++7r05/,7r0)^ - (^27r+7r- + TfO^) (29^7r+5^7r- + a^7r09/,7r0) . 

The next crucial step in the SM is to assume that the physical Higgs boson h^, in addition to 
breaking the electroweak gauge symmetry, also couples to the fermions, generating all the fermion 
masses via the Yukawa interactions, whose strength is given by ratio of the fermion masses mfjr 
over the VEV v, 

yf = V2mf/v, Uf, = V2mf,/v. (4.73) 

Thus the bare fermionic mass-term —nifff — mjiff will be generated from the following gauge- 
invariant Yukawa Lagrangian of the nonlinearly realized SM, 

= -m,(l + ^)^f/(;)/H-m,,(l + ^)^[/(;)/;, + H.c. (4.74) 

The expanded interaction Lagrangian of (4.74) can be directly obtained from Eq. (4.54) by adding 
an overall factor (l + h^/v) . 

The Feynman rules of this nonlinearly realized SM given by Eqs. (4.69) and (4.74) are very 
different from the usual linearly realized SM. But, as they both share the same unitary-gauge 
Lagrangian, they arc actually equivalent and thus both rcnormalizable. In particular, we note that 
the above nonlinearly realized Higgs boson couples to the Goldstone fields with two derivatives, 
similar to the Goldstone self- inter action vertices. This is crucial for maintaining the unitarity of the 
Goldstone boson scatterings vr^ivr"^ — > nvr" and the fermion-anti-fermion scatterings // nvr" 
in the nonlinear SM. For instance, in the latter case, the energy power counting shows that the 
contributions of the EWSB sector in Fig. l(c,d) are of the same 0{E^) as that of the universal 
contact contributions in Fig. 1(a). Hence, the unitarity of the SM will force all 0{E^) terms in 
Fig. l(a,c,d) and Fig. 1(e) to at least cancel down to 0(-E°), i.e.. Fig. 1(a) + (c) -|- (d) + (e) = (D{E^), 
or. 

Fig. 1(a) + (c) + id) = - Fig. 1(e) + 0(^°) = 0{E^) . (4.75) 
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Now we will explicitly analyze the unitarity violation for the scattering //, //' 3ir^ by 
including the EWSB contributions of the SM [cf. Fig. l(c+d)], but without assuming the Yukawa 
interactions of the Higgs boson /i" for fermion mass generation [cf. Fig. 1(e)]. This means that we 
will only include leading contributions from Fig. l(a) + (c) + (d) at 0{E^) [cf. Eq. (4.75)]. Then, we 
will explicitly show how the contributions from a given type of EWSB mechanism such as that in 
the SM [cf. Fig. l(c,d)] could affect the final unitarity bound. We choose the process //, //' — > n7r° 
with n = 3 not only because n = 3 is the simplest nontrivial case where the EWSB effect first 
shows up, but also because the n = 3 case is expected to have the largest EWSB effect on the 
unitarity limit, as explained earlier in this subsection. 




Figure 4: Leading Feynman diagrams of 0{E^) for the scattering //' — > Tr^vr^vr" in the non- 
linearly realized SM. (a) Model-independent universal "contact" contribution; (b) s-channel con- 
tribution due to the Goldstone boson self-interactions from the SM EWSB sector; (c,c') s-channel 
contributions due to Higgs-Goldstone boson interactions from the SM EWSB sector; (d,d') Yukawa- 
type contributions due to the assumption that the same SM Higgs boson is responsible for the 
fermion mass generation. 
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For the n = 3 case, the best unitarity hmit comes from the scattering //' — *■ tt+tt+tt" . The 
relevant diagrams are explicitly shown in Fig. 4. We see that Fig. 4(a) is the model-independent 
universal "contact" contribution (studied in Sec. 4.1); Fig. 4(b) is the s-channel contribution from 
the Goldstone self-interactions in the SM EWSB sector; Fig. 4(c,c') arc the s-channcl contributions 
from Higgs- Goldstone interactions in the SM EWSB sector; Fig. 4(d,d') arc the Yukawa-type con- 
tributions due to the assumption that the same SM Higgs boson is responsible for fermion mass 
generation. Wc then compute the helicity amplitudes for all the individual contributions in Fig. 4, 
at the leading 0{E^) , 



iT [f±f± ^ TT+TT+TT-] (a) 
i7 [f±T± ^ TT+TT+TT-] (6) 
^T [/±f± ^ TT+TT+TT-] (C + C') 
[/±f ± ^ 7T+TT+TT-] id + d') 
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Sl3 



+ 



Sl3 



Sl3 



+ 



where s,- 



(4.76) 

{ki + kjY and s = {ki + k2 + k^Y — ^12 + S23 + S13 > . Summing up all 



contributions in (4.76), we arrive at 



^[/i/'i^ TT+TT+TT-] = {a) + {h) + {c + (^) + {d + d') 



.V2 



S23 



+ 



Sl3 



m 



H 







1 



(4.77) 



which means that the corresponding 2^3 total cross section will behave as 0(l/,s) or smaller. 
According to the unitarity condition (2.26), the scattering f±f'j^ — > tt+tf+tf^ is indeed unitary in 
the SM, as expected. If wc do not assume the Higgs boson to be responsible for the fermion mass 
generation, the diagrams in Fig. 4(d)-|-(d') must be removed, and thus the remaining contributions 
give 

T [/±f± ^ TT+TT+TT-] = (a) + (6) + (c + c') 



= -{d + d'). 



S23 



+ 



Sl3 



S23-m]^ si3-m]j\ 



(4.78) 
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Finally, we evaluate the total cross section numerically using the leading order amplitude in (4.78) 
and plot the results in Fig. 5(a)-(d) in comparison with the imposed unitarity condition (2.26). 

In this figure, we consider four types of initial helicity states, fxf'x G {'t+b+, u^d^, t^T^-, e^V^), 
where the corresponding cross sections are proportional to the fermion masses {ml, m^, rn^, rn^), 
respectively. We have further included the Higgs width in Eq. (4.78) for the actual calculations of 
Fig. 5, but we see that the effects of the Higgs mass and width are significant only in Fig. 5(a) where 
the unitarity violation scale occurs around 3TeV or so; for all other light fermions such effects are 
almost invisible because the relevant unitarity violation scales are substantially above the largest 
values of the Higgs boson mass and width. We see that for all cases in Fig. 5 the inclusion of the 
EWSB contributions makes the cross sections larger for the whole energy range and thus further 
enhances the unitarity constraints, but the corrections to the bounds are small, around 40 — 50% 
or so, i.e., always less than about a factor of 2 for the 2 — > 3 processes. As explained earlier, with 

the increasing number n (> 3) of the final state particles, such corrections will become even smaller, 

1 

approaching unity like [0(2 — 3)]"-! . Therefore, we conclude that the estimated unitarity limits 
based on model- independent, universal "contact" interactions are robust and reliable up to a factor 
of ~2 or much less. 

The Fig. 5(a) has another very interesting feature, namely, the inclusion of the EWSB effect 
pushes the unitarity limit for the scale of the top-quark mass generation down to the region of 

E* = 2.8-3.3TeV, for mj^ = 115 - 800 GeV, (4.79) 

which is below the best bound derived from the scattering tt 2it°' in Table 3. So, we see that 
by including the EWSB contributions (from the minimal SM) , the strongest unitarity limit for the 
scale of the top-quark mass generation occurs at n = 3 , rather than n = 2. Hence, the present 
unitarity limits from 2 — n (n ^ 3) scattering gives better upper limits on the scale of fermion 
mass generation than the classic Appelquist-Chanowitz bounds via 2^2 scattering, for all known 
fermions including the top quark! 

Finally, we note that in addition to Fig. 4 there are three independent 2 — >^ 3 scattering processes, 
// — 7r°7r°7r°, tt+tt^tt^ and //' — tt'^tt^tt'^. We have verified the energy cancellation for each of the 
above amplitudes at 0{E^) when the complete SM diagrams are included. Furthermore, without 
assuming the Yukawa interaction of the Higgs boson hP for fermion mass generation, we compute 
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the amplitudes including the contributions from the EWSB sector of the SM. We find, at 0{E^), 

S23 , Sl3 
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(4.80) 
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Accordingly, we have computed their cross sections numerically and found the unitarity limits for 
these processes to be always weaker than the limits from //' 7r"'"7r"'"7r~ . 

Before concluding this subsection, we would like to remark that the current analysis only makes 

the mild restriction that the longitudinal weak boson V^(= W^, Z£), or its corresponding would-be 

Goldstone boson 7r"(= tt"^, tt^), behaves as a fundamental field below the unitarity violation scale. 

This includes all supersymmetry (SUSY) or non-SUSY models with fundamental Higgs doublets 

for the EWSB. There are dynamical models in which or tt" is composite but the compositeness 

scale is very high, around the GUT scale or above. Examples of this kind of theory include the 

classic minimal top-condensate model [72, 73] and its various viable super symmetric extensions [74], 

and more recently the dynamical EWSB models via a neutrino condensate [75] . In these cases the 

current analysis of the scattering // — nV^ (n7r") will apply. For the case where or tt" is 

composite at a scale of only a few TeV, we may need other types of processes based on effective 

higher dimensional fermionic operators (such as // — {ff)'' with A; ^ 1) for a generic unitarity 

analysis of the scales of light fermion mass generation. 

Recent complementary studies [76] analyzed the low energy phenomenological constraints on the scale of new 
physics via generic four-fermion contact interactions involving the leptonic t — fj, flavor violation. 
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Figure 5: The effect of the SM EWSB sector on the unitarity bound for the scale of fermion mass 
generation, (a). Same hehcity scattering — > tt~^ tt^ 7t~ : the total cross sections including the 
contributions from the SM Higgs sector are shown in solid curves, for rrij^ = 115, 300, 800 GeV, 
respectively. The cross section given by the universal "contact" contribution alone is shown as 
dotted curve, (b). Same as (a), but for the scattering u+d+ — > tt~^tt~^tt~. (c). Same as (a), but 
for the scattering t^u- tt^it~^tt~. (d). Same as (a), but for the scattering e'^u- — > 7r"'"7r"'"7r~ . 
In all cases, including the contribution of EWSB makes the cross section larger and thus tends to 
enhance the unitarity limit, but the effect is always less than a factor of two for the current 2^3 
scattering. Such effects will become even smaller for 2 — > n scattering as n increases above 3 . 
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5. Scale of Mass Generation for Majorana Neutrinos 



Parallel to Sec. 4, this section is devoted to a quantitative analysis of the scale of mass generation 
for Majorana neutrinos. In Sec. 5.1, we will derive the interactions of Majorana neutrinos with 
weak gauge bosons and Goldstone bosons. Then we systematically analyze the scattering processes 
i^i^V]^ ravr" (nV^) and deduce new unitarity limits on the scale of neutrino mass generation. In 
Sec. 5.2, wc discuss the interpretation and implication of our new bounds for the underlying mech- 
anism of neutrino mass generation. Three types of mass mechanisms will be addressed, including 
the seesaw mechanism [20] , the radiative mechanism [21] and the proposal for generating neutrino 
masses from supersymmetry breaking [77] . 



5.1. Majorana Masses and Neutrino-Neutrino Scattering 

Given the observed particle spectrum of the SM, neutrinos can only have the Majorana mass term, 

~ "^"^i/"' ^Tfi^Lj + H.c. , which violates lepton number by two units. With the nonlinearly realized 

Goldstone field U [cf. (2.6)], such a dimension-3 bare mass term can be made gauge-invariant but 

is necessarily nonrenormalizable. As we showed earlier in Eq. (2.10), this uniquely takes the form of 

the Weinberg operator [19] with the nonlinear field defined as $ = -1=^ with $ = [/(O, 1)"^, i.e., 

v2 

£^ = - imt^Lf ^CLf #"'$^'e""'e^^' + H.c. , (5.81) 
2 



„2 



where mj = Cij-^ is the neutrino Majorana mass derived in Eq. (2.11). The neutrino Majorana 
mass term allows us to introduce, for convenience, a (real) Majorana neutrino field xi — X"^) by 
defining 

X = -^{ul + i^l), (5.82) 

where Vj^ = \/2-PlX a-^id = (z^^)'^ = ^/2Prx with Pl,r = ^(1 T 7^) • Using a properly 
normalized notation for the lepton doublet Lj = (^-^Vj^ - , (-Lj^ and the operator (5.81), we arrive 
at the following kinetic and mass terms for the Majorana neutrino field x 1 

^x" = \xj{i^-m,^)Xj, (5.83) 
where m^^ is the jth real mass eigenvalue derived from diagonalizing the symmetric 3x3 mass 
matrix {rriv} . The gauge interactions of Xi which can be deduced from the SM Lagrangian, are 

^? = -^X,7^'X,^°-^V,,^7'^7%W^;7+H.C., (5.84) 
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where = cos 9^] with 0^ the Weinberg angle, and Yij is the Maki-Nakagawa-Sakata-Pontecorvo 
(MNSP) mixing [9] from diagonahzing the Icpton and neutrino mass matrices. In this expression, 
we have used the fact that the Majorana field Xj has vanishing vector coupling. Despite the 
difference between the gauge couplings of the Majorana neutrinos and the massless SM-neutrinos, 
their contributions to the existing physical observables, such as the invisible decay width of etc, 
differ only by a tiny amount of 0(m^/M^) . The real test of the Majorana nature of neutrinos has 
to come from the next generation of neutrinoless double-beta decay experiments. 

Comparing the operator (5.81) with the nonlinear Dirac mass term (2.9), we sec that (5.81) 
also appears as dimension-3 but contains two nonlinear U matrices and violates lepton number by 
two units. Hence, the Feynman rules for the neutrino-Goldstone interactions are different from 
those in Eq. (2.9) though the relevant vertices always have a contact structure with 2 fermions 
and n Goldstone fields vr^i-Tr"^-- • -tt"" (n ^ 1). It is clear that this does not affect our general 
power counting for the amplitude 7[i'j^i']^ nvr"] which takes the same form as 7[ff — ^ n7r"] [cf. 
Eq. (2.17)]. This is why our estimates in Fig. 2 and Table 1 apply to i^ii'j^ — ^ mr°' scattering as well. 

To carry out the quantitative analysis, we deduce from (5.81) the following neutrino-Goldstone 
interaction Lagrangian, in the neutrino mass-eigenbasis,^^ 

= - ^m.^-x,- [A(7r) + i75lB(7r)] X,- , (5.85) 

where A(7r) and B(7r) are given by 

oo 

A(7r)= 3fie(C/22)'-l = J2 ^ 

n (even) =2 



2 ^ 



and 



n-l 

CO -n—l 2 
n(odd)=l i=0 ^ 



2 2^~^ 



(5.86) 



i2k)\in-2k)\ n\ 



even 



n-l 

' n + 1 2^ 



= g (2fc)!(n-2fc-'l)T(2fc + l)(n-2fc) = ^ ' = "^^^ ' ^'"'^^ 

--1 

" g (2A:)!(n-2A:-2)!(2A; + l)(n-2fe-l) " ' = ^ven) , 



*The extension to an initial state like ii'Xj is straightforward. 
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KinCn +i^3nC'li_i 



2n+^-l 




n! 


C n -\- C n 1 
. 2 2 ^. 



with the notation = n\/[£\ {n — 1)1] . Consider the generic scattering process \Xj±Xj±) ~^ 
(7r+)^ (tt-)^ (7r°)"~^^ {£ = 0,1,2,3,---) for three cases: (a) n(even) = 21 ; (b) n(even) > 21 ; (c) 
n(odd) ^ 3 . We compute the hehcity amphtudes 7 Xj±Xj±'i i,i,n — 2£ , from the leading contact 
diagrams [cf. Fig. 1(a)], 



n = even 



7 



_ n n 
Xj±Xj±: ^2 ' "2 ' ^ 



, 1 ITT'uWS 



7[x,±Xj±,t4,n-2£] = ±i"+^e^ 

n = odd: 7[x,±Xj±;iJ,n-2e] = +z»+^e^ "^"^f 
where the dimensionless coefficients 61,2,3 given by 

65' = 2^f.fK^,r., 

6^ = 2'{eif{n-2i)\[Ki,nCi+Ks,nCi_^ 
6^ = 2^{e\f{n-2£y.cLiK2,n. 



(5.88) 



(5.89) 



We further define the normahzed combination for the initial state of Major ana neutrinos with spin-0 
(or spin-1), 

H = ■^[\Xj+Xj+)T\Xj-Xj-)], (5.90) 
corresponding to the case of n =even(odd), and find the corresponding S'-matrix element 



T[in; e,e,n-2e] = ^ ( ^[x,+X,+; A A n - 2£] ^ T[x,_X,-_ -,£,£,71- 2£] ) 



(5.91) 



= V2 7[x^^Xj±;i,^,n-2i]. 

This means that the cross section for the initial state |in) is enhanced by a factor ^/2 . Using this 
we derive the 2 — > inelastic cross section as 

^ ^ (5.92) 



a[m; £,£,n-2£] = 



2Cqj. 



24(n-l)^2n-3 2)! s 

(6;:)2 (e^)2 



(r = 1,2,3). 



Q {d.Y{n-2t)\ ' 

Finally, with the condition (2.26) and = 2\ , we derive the inelastic unitarity limit, y/s ^ E* , 
with 



Et = V 



V 


1 47r 


n-l 









(5.93) 
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where is given by 



'pmax 



rpmax 



rpmax 



23n-47r2n-3 (n - 1)! (n - 2)! 



24("-i)7r2"-3(n-2)! 



n{K2 ^2 



24("-i)7r2"-3(n-2)! ' 
and for n =even, the best Umits are from (3?™^'^, 01^2' 



(n = even, i = ^) , 
(n = even, ^ < §) , 
(n = odd) , 



(5.94) 



Numerically, we find that the bounds 



given by and "^^f^ differ very little, always less than 2%. 
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Figure 6: Scales of mass generation for Majorana neutrinos: the precise unitarity bounds E* from 
Eq. (5.93) are plotted as a function of n (^ 2) , where m^^- = 0.01, 0.05, leV (curves from top to 
bottom). Only integer values of n have a physical meaning. The minimum point of each curve is 
given in Eq. (5.95). 
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With the formulas (5.92)-(5.94), we plot the numerical results of the unitarity limits in Fig. 6 
for three typical values of the neutrino masses mi,j. We find that the best unitarity bound in each 
case to be 

^min = 136, 158, 170TeV, located at n = = 20, 22, 24, (5.95) 

for m^^j = 1, 0.05, 0.01 eV, respectively. Under the same input rrii, = 0.05 eV, this agrees to our 
early generic estimate (Table 1) within about a factor of 1.7. Comparing Eq. (5.95) with Table 3, it 
is striking to note that even though the neutrino mass is about ten million times smaller than the 
electron mass {m^/m^ ~ 10~^), the unitarity violation scale associated with its mass generation is 
only about a factor of 1.3-1.6 higher than that for the electron!^^ For neutrinos such an unitarity 
violation occurs in the multiple nVf^ (n7r") production with the number n = = 20 — 24, while 
for the electron it appears at n = = 11. This leads to an exciting observation that the mass- 
generation mechanisms for both the electron and neutrinos have to reveal themselves at or below 
the scale of ~ 10^ TeV. The scales for the mass generations of other leptons and light quarks are 
much lower, as shown in our Table 2 of Section 4. This not only provides an important guideline for 
model-building, but also indicates where to seek definitive experimental tests, such as the relevant 
astrophyical processes [78] . 

5.2. Unitarity Violation vs. Mechanisms for Majorana Mass Generation 

The unitarity violation scale (2.38) derived from the customary 2^2 scattering indicates that 

the scale of mass generation for Majorana neutrinos could be as high as the GUT scale. At first 

sight, this sounds quite natural if we note that the traditional seesaw mechanism [23] invokes a 

heavy right-handed neutrino whose mass is intrinsically tied to the GUT scale. But, knowing 

there arc other low scale mechanisms for neutrino mass generation [21, 77, 79, 80, 81], we then 

wonder how the unitarity violation limit would just fit into the seesaw mechanism rather than 

other mechanisms especially since the analysis of unitarity violation is universal, independent of a 

particular mechanism. Furthermore, now we have shown in Sec. 5.1 that the unitarity violation scale 

associated with nonzero neutrino mass can be as low as O(IOO) TeV once we open up the inelastic 

multiple gauge boson channels for neutrino-neutrino scattering. What is the physics interpretation 

and implication of our new bounds for the underlying mechanism of neutrino mass generation? 

^^If we only look at the customary limit from 2 — > 2 scattering (cf. Sec. 2.3. 2C and Table 2), we would naively 
expect the mass-generation scale for Majorana neutrinos to be as high as the GUT scale (~ 10^® GeV). The reason 
that the scale is substantially lower will be discussed in the following subsection. 
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We address this issue by analyzing three types of neutrino mass generation mechanisms, including 
the seesaw mechanism [20], the radiative mechanism [21] and neutrino masses induced from SUSY 
breaking [77] . 

5.2.1. Unitarity Violation vs. Seesaw and Radiative Mechanisms 

In our unitarity analysis of Majorana neutrinos (cf. Sec. 5.1), we have strictly followed the same 
logic as for Dirac fermions or weak gauge bosons to define the scale Ajy for generating the Majorana 
mass m^,, i.e., the scale h.^, is the minimal energy at which the hare Majorana mass terra, 

= -^mt^i.i'AL.+H.c., (5.96) 

has to be replaced by a renormalizablc interaction [cf. the general definition given above Eq. (2.12)]. 
A crucial diff"crcncc between the current study and Refs. [44, 66] for the scale of Majorana neutrino 
mass generation is that the neutrinos are not assumed to couple to the same Higgs doublet which 
breaks the electroweak gauge symmetry for their mass generation. We put in the bare mass terms 
by hand for any given type of fermions and derive the corresponding unitarity violation scale as the 
most general and model-independent limit on their mass generation. For instance, if we replace the 
nonlinear field $ in Eq. (2.10) or (5.81) by the usual SM Higgs doublet H, we recover the original 
Weinberg dimension-5 operator [19] which is not a pure mass term even in the unitary gauge and 
implies the assumption that the SM Higgs doublet actually participates in neutrino mass generation. 
To see this more clearly consider, for instance, the traditional seesaw mechanism with a heavy singlet 
Majorana neutrino [20], which gives the renormalizable Lagrangian^^ 

^seesaw = - y'^Lj eH* lyj^ - ^Mru^Cuj^ + R.c. , (5.97) 

and leads to the dimension-5 Weinberg operator after integrating out the heavy Uj^. Because of the 
assumption that the SM Higgs doublet H generates the Dirac mass term —'m^-lJ^^.Vj^ , the resulting 
dimension-5 Weinberg operator invokes the Higgs H for Majorana neutrino mass generation. But, 
by our general proceedure, we put in this Dirac mass term by hand and replace (5.97) by 

^seesaw = " "^Dj " ^^S^i^i? + H-C , (5.98) 

which, after integrating out the heavy z^^, generates the bare Majorana mass term of Eq. (5.96), and 
can be further made gauge-invariant as in Eq. (5.81) with the nonlinear field $. Hence, studying 
the unitarity violation of the scattering Vj^i^i^ mr"' based on a dimension-5 Weinberg operator 
^''By definition, H has hypercharge i. 
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[by integrating out from (5.97)] only probes the scale of the singlet mass Mpc which obeys the 
decouphng theorem [82]^^. But the seesaw neutrino mass is really a combined outcome of both 
bare mass terms in Eq. (5.98) which corresponds to the operator (5.81) after integrating out Vj^. 
Unlike the singlet Majorana mass term for i^^, the non-singlet Dirac mass term in Eq. (5.98) does 
not respect the decoupling theorem [82]. Hence, studying the unitarity violation of the scattering 
v^Vi — > m:"" based on (5.81) probes the scale at which the full seesaw mechanism (5.98) has to set 
in. As will be clear below, even though the singlet neutrino can be still as heavy as AquT; the 
inelastic 2 — > n unitarity requires the seesaw mechanism to start at a scale not much above O(TeV) ! 




(a) (b) 

Figure 7: Majorana neutrino masses radiatively generated at a low scale of O(TeV) ^ Aqut • (a) 
An example of two-loop radiative neutrino mass generation with two extra charged singlet scalars 
(X~, Y ). Here the inserted lepton Dirac masses are given by Eq. (2.9). (b) The Majorana 
neutrino mass contained in the effective mass operator (5.81), after integrating out all heavy scalars 
(X~, Y—) from (a). 

We further note that the most general operator (5.81) need not necessarily arise from a see- 
saw mechanism as in Eq. (5.98). It was demonstrated that even without introducing a heavy v^, 
the desired Majorana mass term (5.96) can be generated from a radiative mechanism at the TeV 
scale [21]. An explicit example of a two- loop mechanism [79] for radiative Majorana neutrino mass 
generation is illustrated in Fig. 7(a). After integrating out the two extra charged singlet scalars 
{X~ , Y ), this generates an effective Majorana mass term for neutrinos in Fig. 7(b) which is just 
contained in the operator (5.81). It is very interesting to note that in Fig. 7(a) this two-loop mech- 
anism also connects the Majorana neutrino masses to the lepton Dirac masses given by Eq. (2.9). 
^^As expected, the resulting unitarity violation scale from the high energy i \\Vi^j^Vi^^) — — > ||Z_LZi, + /i/i) 

scattering [44] is, ~ AtxMr ~ (10^^ - 10^^) GeV ~ Agut, [cf. Eq. (2.36)1, which is 12 - 14 orders of magnitude 

higher than the weak scale. 
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The step from Fig. 7(a) to 7(b) of integrating out heavy singlet scalars {X^ , Y ) is fully parallel 
to the derivation of Majorana neutrino masses from the seesaw Lagrangian (5.98) by integrating 
out the singlet i^^ ! If we consider the special case that the lepton masses are generated by the SM 
Higgs doublet H, we can replace the two nonlinear Goldstone fields <J> in Fig. 7(a) by H and the 
corresponding effective neutrino mass operator in Fig. 7(b) is just the original dimension-5 Weinberg 
operator. Again this is fully parallel to the step of deriving the seesaw neutrino masses from the 
renormalizable Lagrangian (5.97) by integrating out Vj^, which results in the dimension-5 Weinberg 
operator. 

So, once we put in all fermion mass-terms by hand (without making extra model-dependent 
assumptions about their origins), the Majorana neutrino masses generated either from the seesaw 
mechanism or from the radiative mechanism can be formulated as the generic effective mass operator 
(5.81). However, in the seesaw mechanism the new singlet fermion z/^ is known to be as heavy as 
the GUT scale [20] while the new singlet scalars such as {X~ , Y~~) in the radiative mechanism 
can naturally have masses of O(TeV) [21]. Then it is intriguing to ask: how does the unitarity 
violation given by inelastic neutrino scattering I'l^fj^ nvr" [based on (5.81)] constrain the new 
physics scale for neutrino mass generation? We know that this limit must be universal, independent 
of whether the underlying dynamics is a seesaw mechanism, or a radiaMve m,echanism. Starting with 
the general mass operator (5.81), we thus expect that the optimal bound on the unitarity violation 
scale should be neither simply fixed by the mass Mr (~ Aqut) of Uj^ in the seesaw mechanism nor 
by the masses Mx,y (~ TeV) of {X~ , Y ) in the radiative mechanism. This is because neither 
Mr nor Mx,Y are the whole story in each mechanism: there has to be additional Dirac mass terms 
[either the first term in Eq. (5.98) or the lepton mass insertions in Fig. 7(a)] which do not respect 
the usual decoupling theorem [82] . This non-decoupling feature means either a weak-doublet Higgs 
or a nonlinear field U has to be invoked to make the Dirac mass term invariant under the SM 
gauge group, where in the first case the Dirac mass is proportional to the Higgs Yukawa coupling 
(forbidding the decoupling of the corresponding Higgs scalar), and in the second case the nonlinear 
U field links the nonrenormalizable Goldstone-fermion coupling to the Dirac mass. In general, such 
a nondecoupling occurs when we build any fermion mass term involving the left-handed SU{2)l 
doublet fermion together with (i) a weak singlet fermion (giving a Dirac mass), or, (ii) itself (giving 
a Majorana mass). Indeed, it is this common non-decoupling feature that forces the unitarity 
violation scale of Majorana neutrinos to be around O(IOO) TeV, substantially below the GUT scale 
but not much above the scale for electron mass generation despite the fact of m^/m^ ~ 10~^ . 

For the above explicit models of seesaw and radiative neutrino mass generation, this means that 
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our new limits constrain the new physics scale of the leptonic Higgs Yukawa interaction, which must 
invoke extra new fields (such as right-handed neutrinos or Zee-scalars or triplct-Higgs) as needed for 
ensuring renormalizability and generating lepton number violation (although our new bounds do not 
directly constrain the masses of these new fields themselves which obey the decoupling theorem). 

5.2.2. Mass Generation Scale for Majorana Neutrinos in Supersymmetric Theories 

Supersymmetry (SUSY) is attractive for stabilizing the weak scale, but does not explain the 
huge hierarchies in the fermion mass-spectrum or the observed flavor mixings. In fact, it necessarily 
extends the SM with extra three-family superpartners and thus adds new puzzles to the "flavor 
problem" [83]. Furthermore, similar to the SM, the minimal supersymmetric SM (MSSM) [25] also 
does not provide nonzero neutrino masses. Therefore, it is important to extend MSSM appropriately 
to accommodate the massive neutrinos. 

5.2.2A. Scale of Mass Generation for Majorana Neutrinos via MSSM Seesaw 

A minimal extension of the MSSM, for instance, is to introduce a heavy right-handed singlet 
neutrino to form a seesaw mechanism [20] . The MSSM superpotential can thus be generalized 
to 

W = Wmssm + AW^ 

(5.99) 

= [HdLy^e + H^Qy^d - H^Qy^u - nH^Hu] + [-HuLyJi - nMjfi ] , 
where n is a chiral superfield^^ containing a right-handed singlet neutrino 17^ and its scalar- 
partner , while the other chiral superfields have the same meaning as in MSSM. Integrating out 
the heavy n field, we obtain a generic (model-independent) dimension-5 efl'ective operator for the 
MSSM Lagrangian, 

= -^Lf^aLfi7^'i?fe"-'e^^' +H.C., (5.100) 

/ -y„ -I- -I- \ ^ 

where Hu = {^t^ ~ ~j ' = wsin/3, and tan/3 = {Hu)/{Hd) ■ This generates new 

interactions of the Majorana neutrinos with the Goldstone bosons and Higgs particles, 

Ai:r = (uLi^CuLj) {Hu2f + H.C.. (5.101) 



^^The singlet field n can expand in flavor space, e.g., one for each family, and the mass then becomes a 3 x 3 
matrix. In fact, the leptogenesis mechanism [84] for explaining the cosmological baryon asymmetry of the universe 
necessarily requires at least two right-handed neutrinos to ensure the existence of the needed CP-violation phase from 
the seesaw sector (called minimal neutrino seesaw) [85]. 
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This Lagrangian becomes, in the mass eigenbasis, 



m. 



m 



f ^ sin f3 



(5.102) 



where 



vsin/3^° + ^ 



i02 



02 
Pu 



2a 



and 



(5.103) 

(sina, cosa), (s^, c^) = (sin/3, cos /?) and (sg^, S2/3) = (sin2a, sin2/3). The 
fields (vr^, vr^) and [h^, A^, H^, H^) are would-bc Goldstonc bosons and physical Higgs scalars, 
respectively. Here a denotes the usual mixing angle in the neutral Higgs sector of the MSSM [86]. 
As an important feature of the linearly realized MSSM Higgs sector with fundamental Higgs bosons, 
we see that the new interactions in (5.102) contain at most two Goldstonc/Higgs bosons, involving 
the neutral (pseudo-)scalars only. This means that the leading amplitudes of high energy Majorana 
neutrino scattering come from the contact diagrams with a two-body final state, i.e., n = 2 [similar 
to Fig. 1(a)], which are of order E"^ . Hence, the best unitarity limit is expected to be of 0{v^ /m^j) , 
similar to the classic bound derived for the SM case [44]. Also, with (5.102) we immediately deduce 
the scattering amplitudes 



7 
7 



Xj±Xj± 
Xj±Xj± 



w+w- 



H+H- 



+ 0(^0), 



o + o(£;0). 



(5.104) 



whose leading order behavior only starts at ©(E"") , and is due to the exchange of t channel 
charged leptons. Comparing the first amplitudes in (5.104) and (2.37), we note that the reason 
that 7[Xj±Xj± ~^ ^l^l] vanishes at 0{E^) in the present case is due to the cancellation from 
the s-channel Higgs contribution which is absent in (2.37) where we did not assume, a priori, the 
Higgs Yukawa interactions for the neutrino mass generation. In (5.104), the absence of a resid- 
ual 0{E^) contribution to the W^W£ channel becomes evident by inspecting the corresponding 
Goldstone amplitude '^[Xj±Xj± ^ vr+vr^] and the interaction Lagrangian (5.102). 

For the current analysis of the Majorana neutrino scattering into Higgs/ Goldstone boson pairs 
in the MSSM, it is clear that the relevant scattering energy is at the GUT scale. So it is safe to 
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ignore all the masses of the final state Higgs bosons, and accordingly the mixing angle a in the 
neutral Higgs mass matrix. Under this massless limit, we have 0^ ~ /i'^ . Also, we will consider the 
parameter region with moderate to large tan /? ~ 10 — 50 in the discussion below, which implies 
1 and Co ^ 10~^ ^ 1 so that p^l — '^'^ • Thus, for estimating the optimal unitarity 



^ 0.995 



limit, it suffices to set ~ . With these we see that the interactions in (5.102) are greatly 
simplified and reduce to the case of one Higgs doublet, 
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(5.105) 



Prom this we compute the leading amplitudes, 
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(5.106) 
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To derive the optimal unitarity limit, it is desirable to consider the scattering in the spin-singlet 

channel with a normalized final state, \\Xi+Xj+) — IXj-Xj-)] ~^ ^^lZl) + \h^h^)\ ■ Evalu- 

V 2 2 
ating the corresponding s-wave unitarity limit gives, 



< 



m, 



(5.107) 



which is at the GUT scale, in agreement with the bound mentioned below Eq. (2.36). 

It is no surprise that the MSSM unitarity limit (5.107) on the scale of Majorana neutrino mass 
generation from the Weinberg operator (5.100) is as high as the GUT scale. This is because the 
nonrcnormalizability of the dinicnsion-5 operator (5.100) arises solely from integrating out the heavy 
right-handed singlet neutrino I'j^ which obeys the decoupling theorem. The other part of the neutrino 
seesaw sector is the Dirac mass term m^, which is already generated from the renormalizahle Yukawa 
interactions with the fundamental Higgs doublet H^- As a result, the MSSM Higgs doublet does 
participate in the neutrino mass generation, and therefore the minimal scale of the mass generation 
for Majorana neutrinos in the MSSM is already known and is partly set by the Higgs masses of 0(f) , 
which are extremely low and right at the weak scale! This is why including a fundamental Higgs 
doublet in the dimension-5 operator (5.100) causes unitarity violations at the GUT scale which 
could be saturated only by the heavy right-handed singlet neutrino z/^ . In the next subsection we 
will provide an interesting counter example showing that invoking the SUSY breaking sector for the 
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neutrino mass generation would even allow the right-handed neutrino mass Mr to be as low as the 
weak scale and thereby evade the bound (5.107). 



5.2.2B. Unitarity Violation vs. SUSY Breaking Induced Neutrino Masses 

In addition to the SUSY flavor problem mentioned earlier, the MSSM ftirthcr suffers a conceptual 
/i-problem, namely, the supersymmetric /v, parameter in (5.99) is expected to be at the Planck scale 
(Mpi) which contradicts the original motivation of stabilizing the Higgs mass at the weak scale. A 
resolution of this is to invoke a global symmetry G which forbids the vector-like mass term H^Hu 
and is broken only by the SUSY breaking sector. In fact, both the supersymmetric fi and Mr 
terms in Eq. (5.99) are bi-linear and invariant under the SM gauge group, so we wonder why the 
same suppression mechanism would not simultaneously keep both of them around the weak scale. 
This motivates the recent proposal of generating small neutrino masses from SUSY breaking [77] . 
For the current discussion, we will analyze a variation of the model in [77] and effectively generate 
both a dynamical /i-term and a seesaw mechanism at the weak scale. Consider that the SUSY 



breaking comes from a hidden sector at an intermediate scale Ttij ^ ^/vMpi via fields 58 with 
(OS) F = F ftixnfKi vMpi and (05) a = VF ^ vcij . The SUSY breaking is communicated to the SM 
via the gravitational interaction so that the effective theory is cut off at the scale A Mpi . Then, 
take a subset of fields X,Y e ^ , which, together with the field n in the visible sector, are charged 
under G . Thus, the /x-term itself is forbidden by G, but can be dynamically generated from the 
SUSY-breaking operator 

^[X'^HdH^]^, (5.108) 

leading to the dynamical //-parameter, 

f, ^ ^ ^ ^ (5.109) 

^ A A Mpi ^ ' 

Similarly, there are additional operators associated with n, which contribute the following La- 
grangian terms, 

= - '^[Y^nn]D - ^ [YHuLn]p + H.c. (5.110) 

The operators (5.108) and (5.110) can be justified from the i?-symmetry under which {X, Y, n) 
have charges (l, |, |) and L and are i2-singlets. In (5.110) we could have additional operators 
involving the X fields but they are suppressed by higher powers of 1/A. A proper superpotential 
in the hidden sector will ensure that (X, Y) develop VEVs for their F and A components, i.e., 
{X)p = {Y)p = F ^ vMpi « (10^°GeV)2 and {X}^ = (Y)^ = /P. Now, the first term in 
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(5.110) generates a weak-scale Majorana mass for n, k, F/k k, F/Mpi ss v , while its second 
term gives the Dirac mass term ~ v^/F /A ~ v^/vMp[/Mp\ ~ y^w^/Mpi = O(keV) . Hence, a 
non-canonical seesaw mechanim occurs at tree-level, 




A 



A 

F 

T 



diag _ 







(5.111) 



which predicts a right-handed singlet neutrino at the weak-scale in addition to the conventional 

,2^ 



seesaw Majorana masses of ^(^"^J ^^'^ ^^'^ active neutrinos. The light seesaw mass is retained in a 
nontrivial way: the global symmetry G properly suppresses both the Dirac mass and Majorana 
mass Mr simultaneously. Now, we can re-analyze the unitarity violation limit from Majorana 
neutrino scattering XX ~^ ^E^L effective theory with n "integrated out" . Since the effective 

Lagrangian remains the same as Eq. (5.100) in Sec. 5.2.2A, we find the unitarity violation scale to be 
still as high as the GUT scale, despite the fact that the right-handed i/^ mass is actually at the weak 
scale. Indeed, such a SUSY induced neutrino mass mechanism has a "screening" feature that hides 
the real neutrino mass generation scale to well below that revealed from the unitarity violation of 
high energy neutrino-neutrino scattering^^ . However, since the Uj^ is so light, the most effective way 
to probe the scale of neutrino mass generation in this scenario is by the direct production of Vj^ and its 
scalar-partner Uj^, which can be the LSP (lightest supersymmetric particle), and has unconventional 
predictions for cold dark matter^° and for the collider signatures of Higgs bosons [77] . 



^^This "seesaw screening" occurs at tree level rather than loop level, unlike Veltman's screening theorem [87] for 
the SM Higgs sector. 

^°Ref. [88] also analyzed a model similar to [77] for TeV-scale resonant leptogenesis. 
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6. On the Scale of Electroweak Symmetry Breaking 



The weak gauge bosons (= W^, Z^) have the bare mass term, M^W'^^^W" + ^M|Z''Z^, 
which can be made gauge-invariant by the dimension-2 nonhnear operator of the Goldstone kinetic 
term in Eq. (2.8) and is necessarily nonrenormalizable. The unitarity violation in the high energy 
scattering V^^V^^ — > nV£ (n ^ 2) can be conveniently analyzed in the corresponding Goldstone 
boson scattering vr'^^vr"^ — > nvr" according to the equivalence theorem (ET). The power counting 
method in Sec. 2.2 shows that the amplitude for this 2 ^ n scattering behaves as 0( — I , where 



the i?-power dependence is independent of the number of external lines [cf. Eq. (2.17)]. 
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Figure 8: Illustration of the relevant Feynman diagrams for the scattering amplitudes of vr^^vr^^ — > 
n7r" (n ^ 4) which are of 0{E'^/v'^) . (a) The "contact" dia grams; (b) the "t(ii)-channel" type 
diagrams; (c) the "jet-like" diagrams. 



For a quantitative analysis, we derive the pure Goldstone Lagrangian from Eq. (2.8) and arrive 



at 



= ^d>^ie.d^if+ rL (^•^)^'' (^•9^^)'-(^-if)(9^7f -s^vf) 



n-2 



n (even) =4 



n! V 



(6.112) 



where 

If -If = 27r+7r"-h7r°7r° , 



( vf • d^lff -{It -It) (d'^lf ■ d^lt) 



(6.113) 
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{It -It) 



{It -d^itf -{It ■it){d^'it -d^it) . (6.114) 



The cancellation of (vr°)^ terms in the second equation of (6.113) explicitly shows that all vertices of 

type (tt'')" {n ^ 4) in (6.112) are absent, implying that the scattering process Tr^vr" nvr^ {n ^ 2) 
is forbidden at the tree- level. Note that in Eq. (6.112) each interaction vertex contains an even 
number of Goldstone fields^^. In contrast to the 2^2 scattering discussed earlier in Sec. 2.3.2A, 
the scattering Tr^ivr"^ nvr" {n ^ 4) contains additional non-contact Feynman diagrams as shown 
in Fig. 8(b)-(c), which all have the same leading E'^ power dependence. This greatly complicates 
the calculation of the cross sections so that a numerical evaluation becomes necessary. Other 
contributions with internal gauge boson exchanges are at most of 0{g^ /v'^~'^) so that they do not 
threaten unitarity and thus are not shown in Fig. 8. 

To study the 2 — > 4 scattering, the relevant Goldstone interaction Lagrangian is, 
1 _ 1 
QiP 45u^ 

Prom (6.114), we consider all possible 2 — > 4 scatterings, classified into four types according to their 
initial states: 

7r°7r° 7r°7r°7r°7r° (absent at tree level) , 

— >■ TT^ 1T~^ 'K~ , 

— >■ 7r+7r+7r~7r~ (V); 
TT+TT- ^ TrOTrOTrV , 

^ TT+TT- TtV (V), 

^ TT+TT+TT-TT- ; (^-H^) 

tt±ttO ^ tt±ttOttOttO , 

^ tt±tt+tt-ttO (V); 

TT^TT^ — ^ TT^TT^TT^^Tt'^ , 

— >■ TT^TT^TT+TT" (\/)) 

where among each type of processes, the one which has the largest cross section is marked by \J and 
accordingly, shown in Fig. 9(a). In our numerical calculation, each incoming/outgoing Goldstone 
boson field vr" is put on the mass-shell of its corresponding longitudinal weak boson V^. In Fig. 9(a), 
we also depict the corresponding unitarity bound based on Eq. (2.26). We see that among all 
processes the channel tt'^tt^ — > TT"'"TT"'"Tr~TT~ gives the largest cross section and thus the best unitarity 
bound, 

= 3.9 TeV. (6.116) 

In this channel, we observe that the contribution from the contact graph is about a factor of 1/(2— 3) 
of that from the non-contact graphs, and their interference is negative so that the total cross section 

^^This is why in our early estimate of Fig. 2 only the even n values are marked for the unitarity bound on the 
scattering V^^ ^ nV£ . 
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is smaller than the contact contribution by about a factor of 2 — 3 . But, from Fig. 9(b) we see that 
the inclusion of all non-contact graphs only weakens the unitarity bound slightly, from E* = 3.6 TeV 
to E* = 3.9 TeV. This shows that using contact graphs alone can provide a good estimate of the 
unitarity bound even when the non-contact contributions have the same energy-power dependence. 
This is due to the fact that an 0(1) correction to the cross section of the vr^^Tr^^ nvr" scattering 
only affects the unitarity bound by a factor of [0(1)] a^i [cf. Eq.(3.46)], which approaches unity as 
n becomes much larger than 2 . Finally, we may expect to derive a stronger bound by considering 
the isospin singlet channel, / = . For the initial state, the / = combination is given by 
|0i) = Ivr^vr") . For the 47r" final state, using the Clebsch-Gordan algebra we derive the 

a 

following 7 = combination, 

|0f) = [4|7r+7r+7r-7r-) + 4|7r+7r-7rV) + iTrVvrV)] . (6.117) 

V 120 

Thus the scattering amplitude T[0] for the 1 = channel is given by 

7[0] = (Of |T|Oi) = ( 87[+-, + + — ] + 8T[+-, + - 00] + 2T[+-, 0000] ^^^^ 

+4T[00, + + — ] + 4T[00, + - 00] + T[00, 0000] ) , 

where 7[00, 0000] = at tree-level. The corresponding cross section of the 7 = channel is given 

by, for s » > 

^[0] = ^ [ \m\'- (6.119) 

A rough estimate would be to simply take the first five individual amplitudes in Eq. (6.118) to be 
equal, say to the one with the largest value, 7r°7r° — 7r''"7r+7r~7r~. Then, we would have T[0] ~ 
^^T[00, + + — ] ~ 1.37 T[00, + + — ] , which implies 

'^[O] ~ ^ / |1.37T[00,+ + — ]p 

JPS4 (6.120) 

= 1.37^ X 2!^ X c7[00, + + — ] ~ 7.5 o-[00, -t- + — ] , 

where the factor 2!^ is due to the conversion of the phase space integration into the channel 
TT^TT^ 7r+7r+7r^7r~ which contains two pairs of identical particle in the final state. Hence, the 
above estimated enhancement factor 7.5 indicates a reduction of the unitarity bound of 
TT+TT+TT^Tr" channel by factor of (1/7.5) = (1/7.5)8 ~ 0.78 . We see that the estimated unitarity 
bound for the singlet 7 = channel is about 3.9 x 0.78 ~ 3.0 TeV which, as expected, is only 
slightly stronger than that of the best individual channel 7r°7r° —>■ 7r+7r+7r~7r~ shown in Fig. 9(a). 
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Thus, a further precise numerical calculation of the cross section a[0] is expected not to change the 
conclusion. 




Vs (TeV) 

Figure 9: (a) Total cross sections of the scattering processes vr'^vr^ — > vr+vr^^vr^Tr^, ■k~^tt~ — > 
vr'^Tr'^Tr'^'vr", vr^Tr*^ — > tt^tt~^tt~tt^ , and vr^vr^ — > 'it^tt^tt~^tt~ . (b) For vr'^vr'^ — > vr'^vr'^vr^vr", we 
show the individual contributions from the contact graph alone (dashed curve) and from the non- 
contact graphs (dotted curve). The total cross section is depicted by the solid curve. 

In summary, Fig. 9(a), together with the above estimates, shows that the unitarity bound from 
2^4 scattering is about 3.0 — 3.9 TeV, significantly weaker than that from the 2 — > 2 scattering 
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[which is about 1.2 TeV, cf. Eq. (2.33)]. This further supports our earher estimate for nV^ 
(n ^ 2) in Fig. 2 (Sec. 3.2), i.e., tlic strongest unitarity bound for the EWSB scale occurs at n = 2, 
in contrast to the case for the scales of fermion mass generation (cf. Sec. 4-5). 

7. Conclusions 

The requirement from unitarity that the partial wave amplitudes for 2^2 scattering, V^'^V^'^ 
y^'^V^" and //, i^j^i^/^ ^L^^L^ 1 bounded for all energy has imposed nontrivial upper limits 
on the energy scales at which the new physics responsible for the mass-generations of the weak 
gauge bosons {V"' = W^,Z^) and fermions (/, i^^) must become evident. The limit on the mass- 
generation scale for V£ is about 1 TeV [cf. Eq. (1.1)] while the limits for the fermions are proportional 
to /m^^ [cf. Eqs. (1.2) and (2.38)] which is substantially above 1 TeV except for the top-quark 
(cf. Table 2), indicating that the scales of mass generation for all the light SM fermions might 
not be directly accessible by current or foreseeable high energy accelerators. More recently it was 
pointed out [44] that the scattering // nV£ , with n > 2 , could impose much stronger unitarity 
bounds for the light fermions when n becomes sufficiently large; these bounds are so severe that 
they approach the weak scale v = {\/2Gf)~^ — 246 GeV as n becomes large [cf. Eq. (1.3)]. This 
leads to the conclusion that the scattering // — nVl does not reveal an independent new scale for 
fermion mass generation [44] . We further observe that such a bound would unavoidably contradict 
the kinematics [cf. Eq. (1.4)]. We have resolved this contradiction by treating the n-body phase 
space exactly. The present analysis shows that the scattering ff — nV^ (n 2) does reveal 
a separate scale for fermion mass generation. Our bounds from 2 — ^ n scattering with n > 2 
establish a new scale of mass generation for all light fermions (including Majorana neutrinos) and 
are substantially stronger than the conventional 2—^2 limits. 

In the SM, the arbitrary Yukawa interactions of fermions with a single Higgs boson are assumed 
to be the mechanism for the mass generation of all fermions; this is, however, not experimentally 
verified so far. Even if a SM-like light Higgs boson is found to generate the electroweak symmetry 
breaking (EWSB) in forthcoming collider experiments, it could still be completely irrelevant to 
fermion mass generation because such a Higgs scalar is not required to couple to the SM fermions 
via Yukawa interactions by any known fundamental principle. This makes probing the scale of 
fermion mass generation an independent task. The crucial point is to note that all the known 
fermions must join the SM gauge interactions and must also be chiral rather than vector-like under 
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the electroweak gauge group. Therefore, a bare mass term for a SM fermion can be made gauge- 
invariant only in the nonhnear realization of the electroweak gauge symmetry [cf. Eqs. (2.9) and 
(2.10)], so it is necessarily nonrenormalizable and thus will cause unitarity violation in high energy 
fermion-antifermion scattering. The energy where such a unitarity violation happens provides a 
model independent upper bound on the scale at which the new physics responsible for fermion mass 
generation must be manifest. The gauge-invariant, nonlinear formulation of the bare fermion mass 
terms also allows us to efficiently compute the high energy scattering amplitude of // — > nV^ 
from that with the corresponding final-state Goldstone bosons, // nvr", by using the equivalence 
theorem. This is extremely valuable when the final state involves multiple longitudinal weak bosons. 

For quarks and charged leptons, the scale of unitarity violation is shown, as a function of the 
number of final state weak bosons, n, in Fig. 3. These results depend on two opposing factors, one 
of which causes the cross section of 2 — > n scattering to increase as n becomes large, similar to 
Eq. (1.3), and one of which causes the cross section to decrease with n until it eventually overcomes 
the first factor completely. Strikingly, the competition between these two opposing factors makes 
the strongest energy bounds occur at a proper value of n > 2 (except, perhaps, for the top quark). 
The values of n (called n^) which give the strongest upper bounds are summarized in Table 3. These 
new limits depend on the fermion mass and range from about 3.5 TeV (n^ = 2) for top quarks 
to 84 TeV {hs = 10) for up quarks, and from 34 TeV (n^ = 6) for tau leptons to 107 TeV 
{ug = 12) for electrons.^^ Hence, while the masses of the involved fermions decrease by a factor up 
to ~ 3 X 10^ , the bounds weaken by only a factor of ~ 30 in contrast to the customary 2 — 2 limits 
which are weakened by a factor of ~6 x 10^ (cf. Table 2). These unitarity bounds could be affected 
somewhat by model-dependent effects from the EWSB mechanism through self-interactions among 
the y^'s (tt^'s), but such changes are shown to be small, to decrease rapidly with the increase of 
n, and furthermore, to be in the direction of lowering the above limits for a SM or SM-like EWSB 
sector. Such effects indicate that even for the top quark the best bound becomes stronger than the 
customary 2 — 2 limit and occurs at = 3 rather than rig = 2 [cf. Eq. (4.79)]. 

The bounds for the Majorana neutrinos depend on the neutrino mass . With the typi- 
cal values of between 1 eV and 0.01 eV (as suggested by the current experiments on neutrino 
oscillations, neutrinoless double-beta decays and astrophysics observations), the strongest bound 
varies from 136 TeV to 170 TeV with n = rig = 20 — 24. Thus the scale of neutrino mass gen- 
eration is only a factor of 1.3 — 1.6 higher than that for the electron despite the tiny mass ratio 
of m^/m^ ~ 2 X (10~^ — 10"^). These bounds are universal, independent of the details of the 
^^These are the energies which, for instance, could be within the reach of a Very Large Hadron CoUider (VLHC) [89]. 
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Table 4: Summary of the strongest unitarity bound ™rf i'^^ TeV) for each scattering ^^^2 ~^ nVf^ 
(Ci,2 — ^L, f, f, i'l) ^^"i corresponding number of final state particles n = , in comparison 
with the classic 2 — 2 limit -^2^2 (i^^ TeV). For Majorana neutrinos, = 0.05 eV is chosen. 
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mechanism for the neutrino mass generation, even though the underlying dynamics may invoke 
mass parameters of a very different scale such as the right-handed singlet neutrino mass Mr in 
the traditional seesaw mechanism (where Mr ~ Mqut) or in the SUSY-breaking induced seesaw 
mechanism (where Mr ~ O(TeV)), or the charged scalar masses Mx,y in the radiative mechanism 
(where Mx,y ~ O(TeV)). For all these models, our universal new bounds essentially constrain the 
new physics scale of the Higgs-fermion Yukawa interaction, which must invoke extra new fields (such 
as right-handed neutrinos or charged Zee-scalars or triplet-Higgs) as needed for ensuring renormaliz- 
ability and generating lepton number violation (although our new bounds do not directly constrain 
the masses Mr and Mx,y of these new fields as they obey the decoupling theorem [82]). Thus, for 
all known fermions including Majorana neutrinos, the new physics of mass generation must manifest 
itself at a scale below about 170 TeV. 

We see that these bounds are very insensitive to the variation of the fermion masses and when 
the mass value decreases from that of the bottom quark (~ 4 GeV) down to that of the Majorana 
neutrinos (~1 — 0.01 eV), by about 9—11 orders of magnitude, the corresponding unitarity limits 
are only weakened by about a factor of 7 or less. This strong non-decoupling feature for the scale of 
new physics associated with light fermion mass generation is essentially due to the chiral structure 
of the fermion bare mass terms, i.e., the fact that all the left-handed SM fermions are weak-doublets 
but their right-handed chiral partners are weak-singlets or, possibly absent for Majorana neutrinos 
with radiative mass generation. It is this feature that makes the coupling strength of fermions to 
multiple Goldstones (or effectively, multiple longitudinal gauge bosons) proportional to the fermion 
mass, so the decoupling theorem [82] no longer applies. A systematic discussion of this is given in 
Sec. 5.2.1. 
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The summary Table 4 shows that the estabhshed new scales of mass generation for all fermions 
(including Majorana neutrinos) fall in a narrow range, 3 ^ 160 TeV. This not only provides an 
important guideline for model-building, but also motivate and quantify where to look for defini- 
tive experimental tests, such as the proposed VLHC experiment [89] and the relevant astrophysical 
processes [78], etc. 

Finally, for pure gauge boson scattering, Vj^^Vj^^ nVf^ (vr^i-TT^^ nvr") with n ^ 4, the 
quantitative calculations are much more difficult due to nontrivial angular dependences in the 
leading energy terms of the S'-matrix elements, but a comparison of 2 ^ 4 scattering with 2^2 
scattering indicates that the rapidly decreasing phase space factor dominates and the best bound 
remains the customary 2^2 bound of ~ 1 TeV. This picture is further supported by our general 
estimates for the 2 ^ n gauge boson scattering with an arbitrary value of n as shown in Fig. 2 and 
Table 1, where our similar estimates for the upper bounds on the scales of mass generation of the 
quarks, leptons and Majorana neutrinos are consistent with the quantitative calculations within a 
factor 2 or so (cf. Sec. 4-5). Also, the current study of // nV^ (nvr") scattering assumes (vr") 
to be a local field or remain local up to the limit E^j . If (vr") becomes composite much below 
-E^, the pure fermion process // — ^ iff)"' may be useful for a model- independent analysis. 

In conclusion, Table 4 sTimmarizcs all the optimal 2 ^ n unitarity bounds on the scales of mass- 
generation, as compared to the classic 2 — > 2 limits. A short summary of these key results is given 
in [90]. 
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Appendices 



A. Derivation of n-Body Phase Space 

For completeness, we summarize a derivation of the n-body phase space formula (3.43) which is 
used in our analyses in Sec. 3-5. The phase space integral for n particles is 



where P = Pi + P2 ■ For the current purpose, it is safe to ignore the masses of the final state 



particles, i.e., mj /E^ < 1 (j = 1,2, •••,?!). Then we obtain the final formula. 



(7r\ n—1 2 
2) (n-mn-nr 

where s = . 

To derive Eq.(A.122), we carry out the following steps: 

(i) . Replace the delta function in (A. 121) with 

J 5^*\Qn-2 - K-1 - K) 5^*\P - ki kn-2 " Qn-2) , (A. 123) 

and multiply by 

1 = j dSn-2 5{sn-2 - Ql-2) ■ (A.124) 

(ii) . Use the first delta function in (A. 123) to do the integration over kn-i and kn in their center 
of mass frame, i.e., 

d'k„-,d^K, 5(^)(Q,_2-fen-i-fcn) = (A.125) 



/ 



(iii) . Repeat step-(i) by replacing the second delta function in (A.123) by 

j d^Qn-S S^^^ iQn-3 - kn-2 " Qn-2) S^^^ {P - ki kn-Z " Qn-s) (A. 126) 

and multiplying by 

1 = j dSn-3 S{sn-3 - QI-s) ■ (A.127) 

(iv) . Repeat step-(ii) by performing the integration over kn-2 and Qn-2 in their center of mass 
frame using the delta function in (A.124) and the first delta function in (A. 126), 

f^^d^Qn-2 6^'HQn-S -kn-2 -Qn-2) 6{sn-2-Ql-2) = ^ " • (A.128) 

(v) . Continue to repeat the steps (iii) and (iv). 
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The result of steps (i)-(v) is the integral 
J„ = (^)"- Y;,,„£Zil /"rf,,£lZ£i ri,,'JZ£l . . . /'"-",,„_,£2z£Z£=zl , (A.129) 

^2/ 7o S Jo Si Jo S2 Jo Sn-3 

which can be directly evaluated to give (A. 122). 

As an alternate derivation, if we assume the validity of (A. 122) for an arbitrary n 2) , then 
we deduce 

^n+i = ( 77 ) / dsn+i- 



(i)7„ 



' ("-!)!(« -2)! ,^ i3„) 

7r\" s" 



.2/ n!(n-l)! 
which proves Eq. (A. 122) by induction. 



B. Refined Unitarity Conditions 



For completeness we improve the unitarity conditions in Sec. 2.3.1 by keeping the masses of ini- 
tial/final state particles. We first recompute the two-body phase space in Eq. (2.21) without ap- 
proximation and obtain 

327r?7?i 

IJeU^^ZJI = 

'PS2 



(B.131) 



where r]f = r]{s,m^,ml) with (777.3, f^i) the masses of final state particles, and 

I 1 
r]{x, y, z) = - [x^ + + — 2{xy + yz + zx)] ^ . 



(B.132) 



For s ^ mf^mj , r/(s, mf, mj) ~ 1 holds well. In the special case rrzg = 7774 = m, Eq. (B.132) 
gives r/(s, m"^, m?) = [l ~ Am? / s\ ^''^ . For the elastic scattering mi + 7712 + 7774, we have 

^in = ^(*' "^1' "^2) — vi^, fnii ml) = r]f . Thus we can derive a refined form of (2.22), 



r o2 



ce L 



where = ri^^af = y VmW Vin ~ W ~ vi^, "^2) • Accordingly, the unitarity condi- 

tions (2.24) and (2.26) become, 



and 



(B.134) 
(B.135) 
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where the total inelastic cross section is given by <Tinei[2 — > ra] = — / l^ineip ^ . Simi- 

2s?7in Jps„ 

larly we rederive Eq. (2.28) as 



where a^"*^^ = yVinVf^^^ ^^f'^^ ^^id rjf^^ = r](s,m^,ml) . Thus, for 2 — > 2 inelastic scattering, we find 



that the unitarity condition (2.29) becomes 

l<--inel| 



< 



(B.137) 



For the present study in Sec. 3-6, we have s S> (j = 1, 2, • • • , n + 2), so that r/j^, 77^ — > 1 , under 
which all the formulas above reduce back to Sec. 2.3.1. In the traditional case of 2 — >^ 2 scattering, 
a similar effect of r]-^^^ due to the finite masses of the initial/final state was noted before [91]. 
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